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VARIANTS OF GEOMETRIC RSK, GEOMETRIC PNG AND THE 
MULTIPOINT DISTRIBUTION OF THE LOG-GAMMA POLYMER 

VU-LAN NGUYEN AND NIKOS ZYGOURAS 


Abstract. We show that the reformulation of the geometric Robinson-Schensted-Knuth (gRSK) 
correspondence via local moves, introduced in |QSZ14| can be extended to cases where the input 
matrix is replaced by more general polygonal, Young-diagram-like, arrays of the form P. We also 
show that a rearrangement of the sequence of the local moves gives rise to a geometric version of 
the polynuclear growth model (PNG). These reformulations are used to obtain integral formulae 
for the Laplace transform of the joint distribution of the point-to-point partition functions of the 
log-gamma polymer at different space-time points. In the case of two points at equal time N and 
space at distance of order we show formally that the joint law of the partition functions, 

scaled by , converges to the two-point function of the Airy process. 


1. Introduction 

In recent years there have been several important developments in the study of random polymer 
models, interacting particle systems and related statistical mechanics models, which appear to have 
an exactly solvable underlying structure. The object that has been driving these investigations is 
the Kardar-Parisi-Zhang (KPZ) equation and the universality class that this determines. We will 
not embark on an interlude of the ubiquitous nature of the KPZ, but instead we will refer to various 
reviews including |KS92| . |C12| . Even though non gaussian fluctuations, governed by a 1/3 exponent, 
were predicted in |KPZ86] for the systems belonging in this class, the more profound structure begun 
to be recognised much later. Fundamental towards this understanding were the contributions of Baik- 
Deift-Johansson |BD.T99| and Johansson |JQ0] where the combinatorial structure of certain objects 
in the KPZ universality was studied and surprising connections to random matrix theory were 
revealed. In particular, |BDJ99| studied the distribution of the longest increasing subsequence in a 
random permutation and [JOO] studied the distribution of a directed last passage percolation with 
geometric disorder and in both cases the exponent 1/3 for the fluctuations and their convergence 
to the Tracy-Widom GUE distribution was established. 

The two models studied in [BDJ99| . |J00| are fundamentally related to each other via the 
Robinson-Schensted-Knuth (RSK) correspondence. The special (and historically prior) application 
of this correspondence, the Robsinson-Schensted (RS), is a bijective mapping between a permutation 
of integers to a pair of Young tableaux of same shape (one which is standard and one semi-standard), 
in which the length of the first row of each tableau equals the length of the longest increasing subse¬ 
quence in the permutation and more generally the sum of the first k rows equals the maximal length 
of the union of k disjoint increasing subsequences. The RS correspondence was extended by Knuth 
as a bijection between matrices with nonnegative integer entries and pairs of Young tableaux. If 
iwij)i<i<m, i<j<n IS an input matrix, then the length of the first row of the Young tableaux equals 
the passage time in directed last passage percolation given by 

■= max Y Wij, (1.1) 

‘K^Yirn n ' ^ 

(2,j)G'7r 
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where Ilm,n is the set of down-right paths from (1,1) to 

Combined with the works of Okounkov |O01| and Prahofer-Spohn |PS02| . which made connec¬ 
tions to integrable systems and the fermionic formalism, many statistical models of a random growth 
nature where analysed and their asymptotic fluctuations where established, see |J05| for an account. 
However, these developments were restricted to zero temperature models, even though the univer¬ 
sal KPZ/Tracy-Widom fluctuations were also expected for such models in positive temperature. 
In |012| . [C0SZ14| . |OSZ14| it was realised that the geometric lifting of the RSK correspondence 
(gRSK) introduced by Kirillov |Kfll| (who had originally named it “Tropical RSK”-see also [NYOd] 
for a matrix approach) could be used as the stepping stone to make progress in the analysis of fluc¬ 
tuations of statistical models in positive temperature. At the same time the theory of Macdonald 
Processes of Borodin-Corwin |BC14| . based on the theory of Macdonald functions |Mac99| . provided 
another route to this task analogous to the theory of Schur processes of Okounkov-Reshetikhiii 
[OBn3| . 

In this article we work within and further elaborate on the framework of gRSK with the purpose 
of establishing extensions of this correspondence and draw conclusions on the correlation properties 
of a directed polymer model in random medium (DPRM). The DPRM is the positive temperature 
version of last passage percolation (ED and its partition function is given by 

^m,n •— ^ ^ (^■^) 

TT^IIyn.n 

where again H^ „ is the set of down-right paths from (1,1) to (m, n). A distribution of the random 
variables (wij) that leads to an exactly solvable partition (ED is the inverse gamma distribution: 
the variables Wij are considered to be independent and the distribution of each one is given by 
r(ai -|- exp(—l/wij) dwij/wij. A DPRM with such an underlying randomness is now 

known as the log-gamma polymer. It was introduced by Seppalainen in |S12j who, using ideas from 
particle systems, was able to obtain fluctuation bounds of the correct order for log Zn,n- In 
|COSZI4] . |OSZ14| an exact contour integral formula for the Laplace transform of Zn^m was obtained 
via the use of gRSK. This formula was later turned into a Fredholm determinant in [BCR13] from 
which it was deduced that for certain constants ci, C 2 

logZjv.A-CiY (d) 

-atTT^-^ ^gue, 

C2N^'^ N^oo 

where Fgue is the Tracy-Widom GUE distribution characterising the asymptotic fluctuations of 
the largest eigenvalue of a Gaussian Unitary Ensemble of random matrices. 

Besides the above result, which corresponds to the one-point distribution of a polymer starting 
at (1,1) and conditioned to end at (Y, N), one is also interested to know the joint distribution of 
two or more partition functions of polymers starting at the same point (1,1) and ending at different 
locations (fVi, Mi), (fV 2 , M 2 ),... A particular case is when the ending points lie on the same (one 
dimensional) plane {(to, n ): TO-|-n = TV}. Using the space-time coordinate system (a;, t) = {i—j, *+j), 
which is the more natural system for one-dimensional simple random walks, the point-to-point 
partition functions {Zm,n '■ m-\-n = correspond to random walks starting at zero and ending up 
at different points after time N. Thus the joint distribution of these partition functions is known 
as the equal-time-distribution. In the case of last passage percolation the joint law was obtained in 
[■I08| via a discrete version of the polynuclear growth model (PNG) of Prahofer and Spohn |PS02| . 
PNG can be actually viewed as a different construction of the RSK correspondence and its geometric 
lifting will be one of the objectives of the present article. It can also be viewed as an ensemble of 
non intersecting paths, which after appropriate space-time scaling converges to a continuum non 
intersecting line ensemble called the Airy line ensemble . In this paper we will concern ourselves with 
the topmost line in this ensemble, which is known as the Airy process. The Airy process, denoted by 
Ai{-) is a continuous, stationary process whose one dimensional distribution is the Tracy-Widom 
GUE. In [,I08| it was established that when the weights (wy) have a geometric distribution, i.e. 
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P{wij = m) = (1 — q)q^, m > 0 then 


{ '’'N+cpN^/^ x,N -coN^/3 X ~ C1-/V (d) 

C2 


iVl/3 


N—¥oo 


> {Ai{x) — x^} 


xGIR’ 


a^GlR 


for appropriate constants cq, ci, C 2 depending on the parameter q. 


(1.3) 


The analogue result for the joint distribution of the partition functions Zm,n of the log-gamma 
polymer had not been achieved, so far, the reason being that no formula for the joint Laplace 
transform had been obtained. In this work we obtain such formulae. We do this by establishing two 
variants of the geometric RSK correspondence, which we think are of independent interest: 

1. The hrst one is an extension of gRSK in situations where the input array is not a matrix, 
but rather a polygonal shape of the form F. This is done in Section 12.31 by extending the 
reformulation of gRSK in terms of local moves introduced in |OSZ14) . 

2. The second one is the geometric lifting of PNG, which we will refer to as gPNG. gPNG 
appears to put the local moves of |OSZ14| under an intuitive and transparent perspective, 
which may be useful in other situations, as well. The construction of gPNG and its properties 
are described in Section [3] 

Each construction leads to seemingly different integral formulae, cf. Theorems 12.81 and 13.91 for the 
joint Laplace transform of an arbitrary number of log-gamma partition functions corresponding 
to different space-time points. These formulae have a structure that is reminiscent of Givental’s 
integral formula for the GL„(1R)—Whittaker function, but they are nevertheless different (we do not 
attempt to directly check that they indeed coincide). We also point that the formulae are valid for 
arbitrary points, not necessarily at equal time (but some restrictions exist, that exclude the strict 
two-time points, cf. the conditions on (mi,ni) and (TO 2 ,n 2 ) in Theorem 12.121 for progress on the 
two-time distribution in a last percolation model see m and for a physics approach in the case of 
the continuum polymer in [D13b] L Using the Plancherel theorem for GLr!,(lR)-Whittaker functions, 
we are able to rewrite the formula for the Laplace transform of two partition functions at equal 
time as a contour integral, cf. Theorem 12. 121 


To put things into context let us show the formula for the joint Laplace transform of Z(^ 2 ,n 2 )) 

which in the case mi < ni , m 2 > U 2 writes as 


E 




ri' I'M TT n I > 1 TT ’ nj=n 2 -|-l + "^i) 






2 = 1 


^2 ^ T-[‘m2 


/ dfj. s„,(/r) n r(-dj -b M/) TT -4: 


n 


l< 2 <mi 


r(Ai -b fij) 
r(ai -b aj) 


nUmi -1-1 r(Mi + aj) 
(1.4) 


aj 


where £ 5 , £ 5 +-^ are the vertical lines i5 -b iR and A -b 7 -b iR in the complex plane, with A, 7 > 0 such 
that ai < 6 for all i (shortly written as a < A) and aj < d -b 7 for all j (shortly written as d < 5 -b 7 ), 
and s„(A) = Yli<i=^j<n ^(Ai “ Aj)“^ is the so-called Sklyanin measure. Unless otherwise stated we 
will keep the convention that integrals along vertical lines in the complex plane are traced upwards. 


Those familiar with the formula for the Laplace transform obtained for a single partition function 
in [GQSZ14) . [OSZ14| will observe certain similarities to the above formula. In fact, if in the last 
product of gamma functions, J)[r(Ai -b fJ,j)/T{ai -b dj), one replaced the variables with dj and 
combined with the dA integral, he/she would obtain the formula for the Laplace transform of 
Similarly, if one replaced the Ai variables with ai and then combined with the d/r integral, 
he/she would obtain the Laplace transform for .^( 7712 ,^ 2 )- Bearing also in mind that two polymers. 
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one from (1,1) to (mi,ni) and another from (1,1) to (m 2 , 712 ) can only interact (via intersection) 
within the rectangle 1 < * < mi, 1 < j < n 2 (if mi < m 2 and n 2 < ni) , we are lead to see that 
this cross product reflects in a sense the correlation of the two partition functions. 

Having this formula at hand, we work in Section 2] towards obtaining convergence of the joint law 
of two partition functions to the two-point function of the Airy process, in the same spirit as m- 
Our first attempt was to write the above formula as a Fredholm determinant. This, however, is a 
non-trivial task and certainly demands further investigation, in order to understand the (what we 
believe that exists) underlying determinantal structure. We therefore followed an alternative route, 
which consisted in expressing (d as Fredholm-like series. Further nontrivial manipulation of the 
latter expression reveals a block determinantal structure that points towards the Airy process. Even 
though, at this stage we can pass in the asymptotic limit on each individual term of the series and 
identify it with the corresponding terms in the Airy expansion, a cross term arising from the product 
])[ pLj)lT{ai + dij) in (11.411 makes it difficult to justify the uniformity of this limiting procedure. 

We discuss this point a bit more in Section 14.21 For this reason our derivation of the convergence 
of the joint distribution to the two-point function of the Airy process is formal (cf. Conjecture 
01 Section mi. Nevertheless, we believe that our analysis unravels an interesting structure of the 
contour integral formulae. 

Our exact formulae for the joint law of the log-gamma polymer appear to be the first such 
mathematically rigorous derivation for positive temperature polymers. We should, though, mention 
that in the physics literature formulae have been provided for the joint law of the continuum 
directed polymer in relation to the delta-Bose gas with the use of the coordinate Bethe ansatz. 
That was first done by Prohlac and Spohn [PSllaj . [PSllb] and then, again by the use of Bethe 
ansatz but via different manipulations of the formulae, by Dotsenko |D13| . [D14| . Although for the 
same object, it was not obvious that the formulae provided were coinciding. This was later checked 
by Imamura-Sasamoto-Spohn in [ISS13| . Our method is completely different to these approaches, 
since we investigated, instead, the combinatorial structure of the log-gamma polymer. Naturally, 
the formulae we obtain are different but it is worth noting that cross terms consisting of products 
of various gamma functions appear in all situations. Mathematically rigorous formulae for joint 
moments for related particle models, q-Push(T)ASEP’s, have been obtained, see for example |CP15| . 
In principle these can degenerate to formulae for joint moments of polymer partition functions 
for various points (m, n) which would, however, need to share the same coordinate n. Finally, we 
close this review of related results by mentioning the work of Ortmann-Quastel-Remenik |OQR15^ , 
|OQR15b] . There they obtain formulae for the one point distribution (as opposed to our treatment 
of two-point distributions) of the asymmetric exclusion process with flat and half-flat initial data. 
From these, passing to a weakly asymmetric limit they deduce the law for the one-point function 
of the point to plane continuum directed polymer. Their method is based again on the coordinate 
Bethe ansatz. Motivated by related work of Lee [mol, this method allowed them to make a guess for 
the above mentioned formulae. Cross terms (which this time are not given by products of gamma 
functions) appear also in this situation, preventing them from rigorously deriving the convergence 
to Tracy-Widom GOE. 

The article is organised as follows: In Section [51 we start by making a quick review of the 
main objects of RSK (Section 12.111 and of geometric RSK (Section 12.21) focusing on the definition 
and properties of the local moves. We also give a quick reminder of Whittaker functions and the 
properties that we will be using. In Section [5751 we extend the geometric RSK to the case of polygonal 
input arrays and we use its properties to provide the first formulae for the joint Laplace transform of 
several point-to-point partition functions. We also use the Plancherel theory for Whittaker function 
to obtain contour integral formulae for the Laplace transform of partition functions at equal time. 
In Section [3] we construct the geometric PNG and deduce its main properties. In Section [4] we 
elaborate on the formulae for the joint Laplace transform for two partition functions at equal time 
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and we show formally that their joint distribution converges to the two-point function of the Airy 
process. 


2. Geometric RSK on polygonal patterns 

The geometric RSK correspondence was introduced by Kirillov [KOI] as a bijective mapping 

gRSK : (R>o)"^’” ^ (^> 0 )”^"^, 

between positive matrices, via a geometric lifting (i.e. (max, +) replaced by (+, x)) of the Berenstein- 
Kirillov [BK95| implementation of RSK in terms of piecewise linear transformations. Kirillov ob¬ 
served in his paper a number of interesting properties many of which deserve further investigation. 
A matrix algebra framework was set by Noumi and Yamada [NY04| to accommodate Kirillov’s 
theory and this framework was instrumental for |COSZ14] . In |OSZ14) gRSK was formulated via a 
sequence of ‘local moves’ on matrix elements and this formulation was important in exhibiting its 
volume preserving properties. 

2.1. Reminder of RSK, Young tableaux, Gelfand-Tsetlin patterns. TheKnuth 
correspondence takes a matrix {wij)mxn of nonnegative integers and maps it to a pair of Young 
tableaux (P, Q). A Young tableau is a collection of boxes arranged in left-justfied rows of lengths 
Ai > A 2 > • • •. The boxes are filled with the numbers 1, 2,..., n, such that they are weakly increasing 
across rows (from left to right) and strictly increasing along columns (from top to bottom) (this is 
the defining property of a semi-standard tableau). It is useful to think of the input matrix (wij)mxn 
as a sequence of words (I™ii 2 “’i 2 ... 7 ^,“!^)... where denotes a sequence 

of Wij letters j in the word. Each letter in this sequence is inserted in the P tableau via the 
Robinson-Schensted row insertion. That is, a letter j attempts to be inserted in the first row of 
the P tableau by first bypassing all the boxes that contain numbers less or equal to j and then 
occupying the first box that was occupied by a letter greater than j. If all boxes contain letters 
which are less or equal to j, then j is appended at the end of the row. Otherwise, the letter k > j, 
that is replaced by j, is removed from the first row and tries to be inserterted in the second row in 
the same manner. This bumping procedure continues, with a box being eventually added at the end 
of a row in the P tableau, filled with the letter j, and at the same time another box is added at the 
same location in the Q tableau filled with the letter i. One of the consequences of this algorithm 
is that the P and the Q tableaux have the same row lengths Ai > A 2 > • • • and the vector of row 
lengths (Ai, A 2 ,...) is called the shape of the tableau. So P and Q have the same shape. 

Let us note that there are extensions of the RSK correspondence to the case where the input 
matrix (wij)mxn is replaced by general polygonal, Young-diagram-like arrays of the form P. Such 
extensions are related to Fomin growth diagrams |F86) : see also |Kr06| for further discussions around 
these constructions. 

A useful way to record Young tableaux is via the Gelfand-Tsetlin patterns. This is a triangular 
array of numbers {zij)i<i<j<n ^'iid the correspondence to Young tableaux is via 

i 

Zij = tt{of r’s in row j of the Young tableau}, for 1 < j < i < n. 

r=j 

Note that Zni = A^, the length of the i*^ row in the Young tableau. Given the definition of a Young 
tableaux, one obtains the interlacing conditions for the Gelfand-Tsetlin patterns 

Zi+I,j +1 < Zij < Zi+ij, ioT 1 < j < i < n. (2.1) 

The mappings we just described are depicted in the following diagram (for simplicity, we depict 
them in the special case of a square input matrix, i.e. m = n, where both Gelfand-Tsetlin (GT) 
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patterns have a triangular shape) 

( 2 . 2 ) 


^nn ■ • ■ - 2^11 


Zn ... Znl_ 

Even though not depicted, the boxes of Young diagrams (P, Q) are understood to be filled with 
numbers weakly increasing along rows and strictly increasing along columns. In the last mapping 
we glued together the Gelfand-Tsetlin patterns in an u x n matrix. This is possible since the P and 
Q tableaux have the same shape and hence the corresponding patterns have the same bottom row, 
which is placed on the diagonal. The above sequence of mappings explains the reason why RSK can 
be viewed as a mapping between matrices. 

Let us now assume that we have inserted the words (I“ii2’"i2 • • • n“’i") • • • 
and let us denote by {zij) the variables of the Gelfand-Tsetlin pattern after these insertions. We 
next want to insert the word (I“'m>i 2 ’"”*-U 2 ... gg denote by (2)j) the variables of the 

new Gelfand-Tsetlin pattern. After this insertion, the number of I’s in the first row will be 

Zii=Zii+Wml- (2.3) 

This insertion will bump a number of 2’s, which will then be inserted in the second row. The number 
of bumped 2’s will be equal to 

min ( 2:21 - zii,Wmi) = min ( 2:21 - 211 , zn - zn) = - max ( - (221 - ^ii),-(2'ii - ^^ii)), (2.4) 

depending on whether the number of I’s inserted in the previous step (which equals Wmi) bumped 
all the existing 2’s in the first row (which equals Z 21 ~ -Zii) or not. Moreover, a number Wm 2 of 2’s 
will be inserted in the first row, at the end of a sequence of I’s of length max (z 2 i) 2 ’ii) (bumping 
at the same time a number of 3’s) leading to the update 

221 = Wm 2 + max (z 2 i,zii). (2.5) 

The bumped 2’s will be inserted in the second row, thus increasing the number of 2’s there to 

^22 = 2:22 - max ( - (z 2 i - zii),-(zii - zn)). (2.6) 

A similar procedure of inserting and bumping with the remaining letters carries on and can be 
described by analogous set of piecewise linear transformations (local moves). 

2.2. Geometric RSK via local moves and Whittaker functions. The local move 
formulation of gRSK can be viewed as a geometrization (i.e. replace (max, -|-, —) with (-I-, x, /)) of 
the local moves described in the previous section, cf. (E31), (EH), (ESI, EH- Let us recall how it 
works: 

Let an input matrix X = (xij) G (lR>o)"^'". Define In to be the identity map. For 2 < i < n 
and 2 < j < m, define In to be the mapping that replaces the element Xn by Xi-i^iXn and lij to 
be the mapping that replaces the element xij by xij-ixij, while leaving all other elements of X 
unchanged. For each 2 < i < n and 2 < j < m define the mapping lij which replaces the submatrix 

V j 

of X by its image under the map 

(a b\ (bc/{ab + ac) b \ 

\^c dy c d{b-\-c)j' ^ ' ' 


(Wy) 


RSK 


ij Jnxn 


(P) 


211 

222 Z21 


2ll 

Z22 Z21 


(Q) 


■ 2 ^ 711 
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while leaving the other elements unchanged. The upper-left corner of the output sub matrix (EH) 
corresponds to the bumped elements from a row cf. ( 1 ^ . while its lower-right corner corresponds 
to the updated number of elements in a row, cf. (EH- 
For 1 < i < n and 1 < j < m, set 

■= hj o ■ - -o ki- 


For 1 < 7 < u, define 


77™-*+^ 0 - 


i 

< 

m 

0 .. 

-•077™ 

i 

> 

m. 


( 2 . 8 ) 


which corresponds to the (geometric) insertion of the row of the input matrix (or, bearing in 
mind the analogies from the previous section, the insertion of the word). 


We can now give the algebraic definition of the gRSK mapping as 

gRSK := Rn o ■ ■ ■ o Ri. 

Let us illustrate the sequence with an example in the case n = m = 2. Then 

i?l = TtJ = li 2 O ^11 = ^ 12 ) R 2 = 7 r( O 712 = ° ^22 O ^21 = ^22 O ^ 21 - 

gRSK = i?2 O Rl = I 22 O ^21 O ^12- 

and 


gRSK 


/a b\ ^ (a ab\ ^ 121 ^ f a ^ 122 ^ f be/{b + c) 


Let us point out that the mappings Ri defined above can also be written as 

Ri = 92 ° 9 \ 

where 

i j ^1:7— i+1 O ' ' ' O — 1 ° ^ij ^ j 

1 7 + 1.1 O • • • O li_i j_i O lij i > j. 


(2.9) 


ab 

c d) 'Vc d ) ' \ac d j ac ad{b + c) J ' 


( 2 . 10 ) 


Here we used the fact that lij o liij, = liiji o lij whenever \i — i'\ + \ j — j'\ > 2, which is an immediate 
consequence of the fact that, in this case, the mappings lijihiji transform 2 x 2 submatrices which 
don’t overlap. The mappings Qj are also going to be used in the construction of the geometric PNG. 

This representation is related to the Bender-Knuth transformations, cf |K01| . |BK95| : For each 
1 < i < n and 1 < j < to, denote by bij the map on (IR>o)"'^'" which takes a matrix X = (xgr) and 
replaces the entry Xij by 


Xij — (Xi j-i + Xi-I j) I -|- I , 


( 2 . 11 ) 


leaving the other entries unchanged, with the conventions that Xoj = Xio = 0 , Xn+ij = Xi^m+i = 00 
for 1 < 7 < n and 1 < j < to, but Xio + xqi = x~/^i ^ -I- = 1. Denote by Vj the map 

which replaces the entry Xnj by Xn,j+i/Xnj if J < to and l/xnm if j = to, leaving the other entries 
unchanged. For j < to, define 


bn—j+ 1,1 o • • • o bji—i j—i o bjij J ^ n 
^ 1,7 — n-t-l O • • • O byi—i j — i O byij j ^ 71. 


( 2 . 12 ) 


It is straightforward from the definitions to see that the mapping g" : (^> 0 )"^"* e->• ([R>o)"^"* 
satisfies = hj o Vj. 

The next theorem summarises the main properties of gRSK. The first point is due to Kirillov 
[KOI], while points 2., 3. where obtained in |OSZ14| 


Theorem 2.1. Let W = {wij , 1 < j < n, 1 < i < m) and T = {tij ,l<j <n,l<7<TO) = 
gRSK{W). Then 
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fr) 

1. Let Ilm,n be the set of r-tuples of non intersecting down-right paths starting from points 
( 1 , 1 ), ( 1 , 2 ),( 1 , r) and ending at points {m,n — r + 1 ),(m, n), respectively, with r = 
1, m An. Then 


^m—r+l,n—r+l ’ ’ ‘ lm—l,n—llm,n — 


E 


n 


Wi. 




and in particular 


^m,n — E n Wij, 

(*u')G7r 

is the polymer partition function Z„i,n- 

2. Setting tij = 0, whenever {i,j) does not belong to the array ((i,j): l<j<n,l<i< m), 


ti —“t” — 1 

u 


we have 

3. The transformation 

(log Wij , 1 < j < n, 1 < i < m) i-A (log tij, l<j<n,l<i< m), 
has Jacobian equal to zLl. 


The previous theorem was used in |OSZ14j to identify the push-forward measure under RSK 
when the input matrix is assigned the log-gamma distribution. Let us first define the latter 

Definition 2.2 (log-gamma measure). For a sequence of real numbers a = {ai)i>i,a = {ctj)j>i, 
such that ai -|- > 0 for all {i,j), we define the (a, a)-log-gamma measure on W, by 


P(diT) := J| 




r(Q;-i “1“ ^ 


(2.13) 


Introducing the notation (we will also be using A for minimum and V for maximum) 

TzAn—1 . 


n 


m/\j-l , 

r—0 ^m—r,j — r 


n 


r—0 


and Tin := 


^m—r,j — r—l 


n i/ \ 11 — j. 1 

r=0 ^i—r,n—r 

n (i-l)An-l , 

r—0 ^i—r—l,n—r 


for j = 1, ...,n and i = 1, we have that, cf. |OS7d4| . Corollary 3.3. 

n m 

p o gRSK- ‘ (dT)=n (n C')(n u„“‘) n 

ij \ j/ j—i 


dt 


u 
t - ’ 

b?,'} 


(2.14) 


where in the above products the indices i, j run over i = 1,..., m and j = 1, We call the vectors 
(TTnj)j=i,...,m (Tin)i=i....,Tn the (geometric) type of the (geometric) P and Q tableaux, respectively. 
This definition is the geometric version of the type of a Young tableau, which in terms of Gelfand- 
Tsetlin variables is defined to be the vector ( ^ records the number 

of letters 1 , 2 ,... in the tableau. 

Let us now assume, for the sake of exposition, that the input matrix IT is a square matrix, i.e. 
m = n. Integrating the measure P o gRSK“^(dT) over the variables tij, i ^ j, we obtain the law of 
the bottom row of the Gelfand-Tsetlin patterns cf. (12.21) (which encode the geometric lifting of the 
shape of the corresponding tableaux. This is given (in terms of Whittaker functions) as follows 

n , , 
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. «11 

/\ 

• Z22 • Z21 

/\/\ 

■ Z33 • Z32 • Z31 

/\/\/\ 

• Z44 • Z43 • Z42 • Z44 


Figure 1. This pictures shows a Gelfand-Tsetlin pattern and the arrows depict di- 
agrammatically the potential in the integral form for the Whittaker function, which is 
“ Z,alZh- 


Whittaker functions have an integral formula due to Givental 


xn(") 

ai,. 


, (^15 ■ • ■ 5 ^n) — 




n 


Zi\ ■ 






<j<i<n 


n 


dz,;. 


l<j<i<n—1 


where Z{x) = {(zij)i<j<i<n G : Znj = Xj, forj = If we change variables 

and define ip^i,---,a„ix) = where Xi = log Zi for i = then the functions 

4’^!,...,Ctrl(x) are eigenfunctions of the open quantum Toda chain with Hamiltonian given by 




n— 1 


^Xi+i-Xi 


with eigenvalue — Moreover, they play an important role in the theory of automorphic 

forms |G06| . In our setting they appear more naturally as generalisation of Schur functions having 
the form of generating functions of geometric Gelfand-Tsetlin pattern. The latter are to be thought 
as triangular arrays in the spirit of the standard Gelfand-Tsetlin pattern, with the difference that 
the interlacing conditions m are violated. Instead, the hard boundary conditions of relations m 
are ‘softened’ via the potential 


- S 








which penalises exponentially patterns that violate the interlacing conditions, see Figure [T] 

One of the reasons that the log-gamma measure is amenable to analysis is the fact that Whittaker 
functions have tractable Fourier analysis. In particular, there exists a Plancherel theorem [STS94| . 

[KLOT]: 


Theorem 2.3. The integral transform 



i=l 


dxt 

Xi 


defines an isometry from T2((R>o)”i nr=i T 2 *^™((tlR)", s„(A)dA), where is the 

space of L 2 functions which are symmetric in their variables (the variables Ai, A 2 ,... are unordered), 
i = if—l and 


Sn(A) 


1 

(27rt)"n! 


nr(A.-A,)-\ 

i¥=3 


is the Sklyanin measure. That is, for any two functions f,g G L2((R>o)”i 0"=! holds 

that 


'(R>o)’ 


f{x)g{x) 


dxi 


i=i 



/(A) 5 (A)s„(A)dA. 
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Naturally, this Plancherel theorem will play a crucial role in this work, as well. We will also use 
an integral identity involving Whittaker functions with importance in the theory of automorphic 
forms (see [B84]) 

Theorem 2.4 (Stade). Suppose r > 0 and A, u G C", where 5R(Ai + i/j) > 0 for all i and j. Then 

f f[^= + A,). 

J(R>o)" „■ 

This identity was conjectured by Bump in |B84| and proved by Stade in |St02| . An alternative 
proof using the gRSK was provided in |OSZ14] . Corollary 3.7, which also showed that Stade’s 
identity is the analogue of the Cauchy identity for Schur functions. 

When using the Plancherel-Whittaker theorem we will have to check that certain functions, which 
are expressed as products of gamma functions, are integrable. For this will use the following 
asymptotics of the gamma function, cf [OLBCIO] . formula 5.11.9 : 

lim |r(a-ht6)|e5l^l|5|^-“ = \/^, a > 0, 6 G R. (2.15) 

b—¥oo 

2.3. Extension of gRSK to polygonal input arrays. In this section we will extend 
the gRSK to the case where the input array is a polygonal array. Young-diagram-like of the form 
P, rather than a matrix. Let us start by considering pairs of positive integers {mi, n^)^=i, 2 ,...,fc with 
mi > mi-i and > ni and rectangular arrays = (wij : mi-i < i < mi,l < j < ni). 

By convention, we have that mo = uq = 0. The general input pattern will consist of an array 
W = lJf=i created by stacking the rectangular pattern below (see Figure [5] for an 

illustration). We define the index set of such a polygonal array W as 

k 

ind{W) := |J {(i, j): mi-i <i<mi,l<j < ni) (2-16) 

1=1 

We now come to 

Definition 2.5 (gRSK on polygonal arrays). We define the geometric RSK on the polygonal 
array W = \_^i^i inductively as follows: 

(1) Set := gRSK{W^^'') G is then the rectangular output provided by gRSK, as 

this is described in Section \2.‘A 

(2) Create the array |J by appending matrix at the bottom of the output matrix 
obtained in the previous step. We then define ;= Rm^ ° Rm^-i o • • • o R^^+i (W(2) IJT(i)). 
The crucial point, here, is that the transformations Ra, mi -I- 1 < a < m 2 , act only on the 
entries {i,j) of |J T^^'^ for which i — j > mi -|- 1 — ri 2 - That is, they only transform the 
entries which lie on or below the north-west diagonal starting from entry {mi 1 , 712 ), while 
leaving all other entries unchanged, see Figure\^ 

(3) We can, now, define inductively T^^'> := R^t o Rme-i o • • • o {W^ Ur=i for i < k 

and, finally, gRSKfW) := T^^fi 

The next proposition summarises the fundamental bijective properties of gRSK on polygonal 
arrays and it is the analogue of Theorem 12.11 

Proposition 2.6. Consider a polygonal array W = ut with = (Wy , 1 < j < 

ni, mi-i < i < mi) and T = {tij, {i,j) G ind{W)) = gRSK{W), where ind{W) is the index 
set of the polygonal array W defined in (12.161) . Then 

(r) 

1 . Let Ilrn{n be the set of r-tuples of non intersecting down-right paths starting from points 
(1,1), (1, 2),..., (1, r) and ending at points {m, n — r ..., (m, n), respectively, with (m, n) G 
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’ (6, 6) 













’(8,4) 






' (11,3) 


- 4 


(4,10) 


Figure 2. This figure shows a rectangular array W = [J (J VF® [J 
which is a concatenation of four rectangular arrays. Their bottom right corners are marked 
by points (me, ne)e=i ,...,4 = (4,10), (6, 6), (8, 4), (11, 3) 


































(6,6) 








. fa /l^ 







, m 




(4,10) 


Figure 3. The red contour shows the part of the polygonal array that will be 
transformed after application of the gRSK on the two top rectangular arrays, i.e. 
gRSK(VF*-^^ IJ W^^). The blue dots show the entries that will be modified by the insertion 
of the seventh row (ui( 7 ^i), u'( 7 _ 2 )W( 7 , 3 ), W( 7 , 4 )) via application of transformation R 7 . From 
this it is clear that R 7 does not see the full polygonal array and it would lead to the same 
output, as if only the subarray , 1 < i < 7, 1 < j < 4) existed. 


{(me, ne ): £ = 1, ..., fc} and r = m An (see Figure{^. Then 


tm.— 


m—r-\-l,n—r-\-l * * 'l,n—l^m,?! 


im,— 1 .n— 1 — 


E 


n 


vtj , 


{7n,...,7rr)£ni^]n J ^ U *' * UtT^- 


and in particular 


tn 


E n 


Wij, for £ = l,...,k, 

are the polymer partition functions Zmt,ni ■ 

2. Setting tij = 0, whenever (i,j) does not belong in the index set of the polygonal array, we have 

ti — 4,j ti^j — 1 


{i,j)^ind{W) ^ {i,j)£ind{W) 


tij 
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(6, 6) 













'(8,4) 










(4,10) 


Figure 4. This figure shows three non intersecting, down-right paths in The 

partition function corresponding to this ensemble equals t 6 , 6 l 5 , 5 l 4 , 4 , where (Uj) = 

gRSK(Uti 


3. The transformation 

(logiUy , (i, j) e ind{W)) {logUj , {i,j) e ind{W)), 
has Jacobian equal to zkl. 

Proof. Part 3. follows immediately from the volume preserving properties of the local moves. The 
proof of part 2. has minor modifications of the proof of part 3. of Theorem 12.11 Since a similar 
property holds for the geometric PNG, whose proof will require slightly more modifications, we 
prefer to defer the details till Proposition [231 

The proof of Part 1. proceeds by induction in k. For k = 1 this is the same statement as that of 
Theorem mi Part 1. Assume, now, that the statement is valid for fc — 1, that is, 

k—1 

:= gRSK( y IPW) 

^ e=i 

satisfies the statement of Part 1. Let 

T = r« := o o... o 

Observe that a transformation Ra acts only on the entries {i,j) with i — j > a — jf{a), where 
j*(a) = max{j : (a, j) G ind{W)}. This means that when Ra inserts a line {waj, 1 < j < j*(o))) then 
it will act as if there was only a rectangular pattern {wij, 1 < a, 1 < j < j*(a)), ignoring the rest of 
the polygonal pattern. This, in particular, holds for a = mk-i + 1, ...,mfe. Moreover, the insertions 
of (wij) with j > rife (that have taken place before the insertions Ra for a = mk-i + 1 , ...,mfe) will 
not have an effect on the entries {tij : i — j > ruk-i — rife) of T (see also Figure |21). Both of these 
statements can be seen as a consequence of the Bender-Knuth transformations (12.1111 . (I2.12|l . which 
imply that when an entry Wij is inserted, it leads only to updates of the entries of the array along 
the north-west diagonal starting from (i, j). 

The claim now follows by combining the induction hypothesis with the statement of Theorem 
nn Part 1, applied to the rectangular array {wij, 1 < r < rrife, 1 < j < rife). □ 

It will be useful, for later purposes, to state separately a proposition, which defines and identifies 
the type, cf. (12.141) . of general polygonal patterns via the gRSK: 

Proposition 2.7. Let a polygonal array W = Ufci with = (wij , 1 < j < Ui, m^_i < 
i < me) and T = gRSK{W). We have that 
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• For every j, we set if{j) := inax{i: (i,j) G ind(W)}. Then 

(2.17) 


(2.18) 

Proof. We will prove (12.171) . the proof of (12.181) being similar. When k = 1 the polygonal array 
is a rectangular array and the claim follows from the properties of the standard gRSK, cf. point 
1. of Theorem O Assume, now, that k = 2. Since after the insertion of the entries (ty) 

of gRSK(W^^^ IJ are identical to the corresponding ones of gRSK(lR^^^), claim (12.171) hold 

immediately for ^2 < j < ni. For 1 < j < n 2 the entries (ty ) involved in the right hand side of 
(I2.17|) have indices satisfying i — j> m 2 — n 2 but as already mentioned the values of these entries do 
not depend on (wij) with j > n 2 , which means that these entries would be identical to the entries 
of gRSK((wy : 1 < i < m 2 ,1 < j < ^ 2 ))- Bearing this in mind, the validity of (12.171) for 1 < j < n 2 
is a consequence of the properties of standard gRSK (Theorem 12.II point 1.). The same argument 
continues when fc > 3. □ 


n 

2=1 


W.i = Ti 




n ** \j 
r=0 




n;*i 




0 




For every i, we set jf{i) := max{j: (i,j) G ind{W)}. Then 

pri.(i)Ai-l , 
llr-=0 


F (») 

rCy = Li 

1=1 


TTl*(dA(j-l)-l J. 
llr=0 


2.4. Multipoint distribution of the log-gamma polymer. We are now ready to 
use the formulation of gRSK on polygonal arrays, in order to obtain expressions for the multipoint 
distribution of the log-gamma polymer. The first step in this direction is to determine the push 
forward measure under gRSK of the log-gamma measure on a polygonal array W = Ufci 
with = {wij : mi-i <i < me,l < j < ne) with {mi,ni)i=i^ 2 ,...,k such that mi > mi-i and 
ni-\ > ni- By convention mo = ng = 0. To do this, we first recall the definition of an (a, Q;)-log- 
gamma measure on such a polygonal array W: 


P(dW) := 

(i,j)Gind(W) 

Rewrite the (a, Q;)-log-gamma measure as 


1 (Oi + aj) L 


1 


m *.( 1 ) 

■=(*>')= n r(o+d)iivii 

(iJ)G ind(W) ^ ^ J = 1 i=l 


F(i) 


n(n»«)’'n(n 


w. 


X exp 


(- E n 


i=i i=i 

dw- 




(■i ind{W) 


ind{W) 


W. 


^3 


By Proposition 12.61 and Proposition 12.71 and by changing variable {wij: {i,j) G md(lV)) 1 —>■ 
(tij: {i,j) G md(lV)) = gRSK(lT) we have that the push forward measure is 


Po(gRSK)-i(dr) 


n 

(i,j)G ind(W) 


1 

r(ai -I- dj) 


(r.(1)j) n (r.I.I*)) 

j=l i=l 


xexp(-£(T)) n 

(i,j)G ind(W) 


(2.19) 


We now arrive to our first formula for the joint law of the partition functions: 

Theorem 2.8 (k-point distribution). Let (mi, ni),..., (wfc, rife) be k pairs of positive integers 

such that mi > m^_i andni-i > ni and consider the point-to-point partition functions Z(^rni,ni), * = 
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1, k of an (a, a)-log-gamma polymer. The Laplace transform of the joint distribution is given by 


— u-\Z(. 


-UkZ(^„ 


'^fc) 


I 


— Uitrr 




Po {gRSK)-\dT), 


where P o [gRSK) ^ is the push forward measure under gRSK of an (a, a)-log-gamma measure on 
the polygonal array W := ni) (mknk)) determined by the corners (toi, ni),(m/c, n/c), see 

mm . 

Proof. The proof follows immediately from the fact that if T = (tij) = gRSK(lT(mi,ni),...,(mfe,rafc))) 
then by Part 1. of Proposition 12.61 


^{rnk.rik))^ 

and the identification of the push forward measure P o (gRSK)“^(dT) as in (I2.19p . □ 

We will now use the Plancherel Theorem l2.3l for Whittaker functions to rewrite the above integral 
as a contour integral in the case of fc = 2 . We first introduce some notation 

Definition 2.9. For a sequence of complex numbers a = (oi, 02 ,...) and integers p, r we will denote 
the function 


^p,rM ■= n r(it; + aj), forwGC. 

p+l<z<r 


When p = 0 we will simplify the notation by denoting the corresponding function by Ff(w). 

We will also make use of the following theorem proved in |COSZ14j . [OSZ14) . which expresses 
the Laplace transform of a single partition function as a contour integral 

Theorem 2.10. Let Z(^jn,n) be the partition function from (1,1) to with m > n, of an 

(a, a)-log-gamma polymer such that d < 0 and a > 0. Its Laplace transform is given in terms of a 
contour integral as 


— uZ, 


(m,n) 




d/4 Snip) n n ~ 


l^j 






U “4 


Remark 2.11. If the parameters a, a are such that + ^) > 0 and 
i5 G R, then the above formula writes as 


FZ(Pj) 

Fm(di) 

-d + (5) >0, for some 


E 




n 

d/4 Snip) n n 

i=i 




where (.g is the contour along which the real part is constant and equal to 6. 

We can now state 

Theorem 2.12 (2-point contour integral). Let Z(^n%i,nk)^ Z(^n% 2 ,n 2 ) bs the partition functions from 
(1,1) to (mi,ni) and (m 2 ,U 2 ), respectively, of an ia,d)-log-gamma polymer. Assume that mi < m 2 
and ni > n 2 . Moreover, without loss of generality, assume that mi < ni. Let 7/2 > 5 > 0 arbitrary, 
such that the parameters (ai)i=i, 2 ... OLnd idj)j=i^ 2 ... satisfy |a| < S, |d — 7 I < S, and let be 

vertical lines in the complex plane with real part equal to 5 and (5 + 7 , respectively. The joint Laplace 
transform is given in terms of a contour integral as 
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in the case that m 2 > n 2 , then 


3 ^ 2 ^m 2 ,T) 




™i c& /"xA 

dAs^,(A) n r(-a, + A. 0 n;;^ 


i=i F4,„^(a*) 




^2 Mj po: / A 

d^^sM n r(-d,+M/)n ' 

i<ij '<"2 


n 

l<z<mi 


r(Ai -i- fij) 

r(ai -I- aj) ’ 


= 1^2 

( 2 . 20 ) 


• in the case that m 2 < n 2 , then 


^ — UiZmi ,711 ”'^2-^77 


/ mi 

^ dAs™,(A) n r(-a, + A,0n 

^ l<i,i' <mi i—1 


{U1IU2) ^^fZ.ni{>^i) 
’ (wi/u 2 )-“‘ F“^ 


'(V)" 


d/iSm2(A*) n r(-ai-|-/ii/) 




mi + l<i<?7i2 

l< 2 ^<m 2 




“ 2 “' F 4 (ai) 


X r(—Ai-i-^i'), 

l<z<mi 

l< 2 ^<m 2 

where the contours is, ^ 5 ' have real parts 6,6', respectively, with 6' > 6. 


( 2 . 21 ) 


Proof. We first treat the case m 2 > n 2 and for simplicity we assume that m 2 > n 2 (the case 
m 2 = n 2 requires only minor notational modifications). From Theorem 12.81 we can write 


g UiZrni ,ni "^2^7712,712 

1 _ / 

-I- O'™ ! I 

i=i 


^ / g-«lt™i,ni-«2t™2,n2 P O (gRSK) ” ^ (dT) 


n*jF(ai + dj) 


ni m 2 

IT (^**0)j ) II 


X exp ( - ^ 


2=1 

ii —— 1 


t 


e til ^2^7712,' 

dti 


n 


^ij 


' ii,j)£ind{W) 


We integrate first over all variables (ty) except {tmi,ni,tmi-i,ni-i, and also do the 

integrations over the variables (tij : i — j < mi — ni) separately from the integration over the 
variables (tij : i — j > mi — ni). We can then rewrite the above integral as 


mi —1 


nr(a. 


0.6 


qiR>o)”*'- 




r—0 


d^mi —r,ni —r 

~t ’ 

^mi—r,ni —r 


( 2 . 22 ) 


where denotes the vector (tm-i^m,trm-i,m-i,—ti,m-rm+i), is a GLmd^)- 

Whittaker function corresponding to the triangle (tij: i — j < mi — ni), i.e. (recall also notation 
Ti j^(^i') from Proposition 12.71) 


vi/(”n)(' 7 '\ ) 

a \ mi ,ni J 

r —O' 

= /n(^Li.w) exp( 

i—1 


— l.j H” — 1 

{i,j)Gind{W) 

i—j<.mi—ni 



E 

ind{W) 
i—j—mi —ni 


tij / 


n 

(z,j)G ind{W) 
i—j<.mi—ni 


(2.23) 

^tij 



















16 


V.L. NGUYEN AND N.ZYGOURAS 


and is a Whittaker-type function corresponding to the pattern {tij '■ i — j > mi — ni), 

defined by 


/ rii m 2 

rr n 

j — l 


i: 


{i,j)G ind{W) 

—ni 


—1 
tij 


E 


{iJ)G ind{W) 
i—j—mi —ni 


tin 


(2.24) 

{i,j)Gind(W) 

2 —j>mi—ni 


We will now use the Plancherel Theorem l2.3l for Whittaker functions to rewrite the integral in (j2.22p 
as a contour integral. First, using the Plancherel theorem, we write (j2.22p as 



X 


mi — 1 




:,-Uitmi,ni 




r—0 


—r,ni —r 
tm.i —r.ni —r 


r,ni—r 


mi —1 


„ — U2trt 


^ Q;,a v^mi ,ni / ^A V‘'mi,ni/ J_ 


r—0 


tmi — r,ni — r 


(2.25) 


Moreover, we shift the contour of integration from (ilR)"*i to (€ 5 )™^ We do this to avoid the 
encounter of poles later on. In order to justify both (12.2511 and the shift of contours, we need to 
check that the two factors inside the parentheses belong to ; 5^1 (A)dA). This will be easier 

to check once we derive some more explicit formulae for these factors. The first integral in (12.251) 
can be computed by Stade’s identity 12.41 and equals 

Y[ r{-a, + Xe). (2.26) 


Notice that the above quantity belongs to L^(( 6 lR)™i; Smi (A)dA) , since by the asymptotics of 
Gamma function (12.151) . we have that for all values of Ai,..., Ami with large imaginary part 


u(A) 




|r(-ai-b Ai/)p Ri exp I ^ ^ |Ai - Ai'l - TTTOi ^ |Ai 


l<2,z'<mi 


<mi 
mi \ 


< exp -TT^ |A,| , 


(2.27) 


which decays exponentially. 

If we, now, expand the functions and „^) in terms of their integral rep¬ 

resentations (I2.23|) and (I2.24|) . we recognise that the second integral in (12.251) . is the Laplace trans¬ 
form of the point-to-point partition function ^ (Ai,..., A^i, 0 ^ 1 + 1 , ■•■am 2 ) di,..., q:„j)- 

log-gamma polymer multiplied by the normalising factor 


r(Ai-|-Q:j) T{ai+aj). 

l<2<mi mi+l<i<m2 


We can, therefore, write 


/ mi-l , 

.^_n ^mi-r,ni-r 


= n ua.+a) n 


r—0 

r(ai -b aj) E 


-u2z2:' 


l< 2 <mi 


mi+l< 2 <m 2 


(2.28) 


To make this formula more transparent, compare with (12.221) . set ui = 0 and recall that tm 2 ,n 2 is 
the partition function Zm 2 ,n 2 - 
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When 7712 > 712 we may use the formula provided by Theorem 12.101 applied to the Laplace 
transform of ^ < 5i((5 + 7 — d) and 5R((5 + 7 + a) > 0, in order to write (12.2811 as 

n r(A,+o,) n r(ai + dj) / dfj.Sn2 (m) n r(-d/ +fij) 

mi + l<i<m2 (^* 5 + 7 ) ^ "^^2 

-fr UT=\ r(M,- + A.) UTJm,+i + a.) 

"" M UTJi r(«i + A,) r(d, + «,) 

(2.29) 

Finally we insert (12.291) and (12.261) into (12.251) we arrive to formula (I2.20p . In a similar way we 
can arrive to (|2.21|) . in the case m 2 < n 2 - The only difference would be that we will need to adapt 
Theorem l2.1()l to apply it when m < n, but this is straightforward by transposition of the rectangular 
array. 

It only remains to check that (12.281) or equivalently (12.291) belongs to L^((7.IR)'"7 Smi (A)dA). We 
prove this separately, in the following lemma, as it requires a particular combinatorial analysis. □ 

Lemma 2.13. The quantity (|2.28|1 (or equivalently (12.2911 ) belongs to 

Proof. We work on the integrability of (I2.29P and start by using the Cauchy-Schwarz inequality 
(bearing in mind the easily verified fact that the Sklyanin measure Sn 2 {h) is ^ positive measure): 



" 2 -l-l<j<rai 



dfj.Sn2{h) n 


< 


n2 mi m 2 

xnu-^^^-“^')nr(M7+A.) n r(/7,+a.) 

j—l i—l i—mi-j-l 


[ dXSmiW n 

|r(A*+dj)|^/ dfiSn2{h) n 

+ Ai) 

' n 

2 

T{qij + ai) 




mi + l<i<m2 



n2 + l<j<ni 


'(is+^Y‘ 


dil Sn2{h) 


- I2l=i(Uj-aj] 


n 

l<i:f '<"2 


T{-af +fij) 


The second integral with respect to dfj, is independent of the A variables and finite for the same 
reason as in (12.271) . To check the finiteness of the rest of the integral we use the asymptotics of the 
gamma function (12.151) . Since the variables (Ai),(/ij) have a constant real part and since we will 
only be using the exponential part of the asymptotics of the gamma function (12.151) . which only 
involves the imaginary part, we may assume that the real part is zero. We then have 


^mi (A) Sn2 ih) 


n ir(A.+d,)i" n 


'n2-\-l<j<ni 


l< 2 <mi 


T{fij + Aj) 


2 

mi + l< 2 <m 2 



r{fij + ai) 


2 


|A*-Ai/| + | 1^*1 

n2+^<j<ni 

-TT |A,+/ 7 j|- 7 r Y IYjI)- 

mi + l<z<m 2 

l<_7<n2 l^i^'^2 


Since the last expression involves only absolute values, we may further assume that A’s and /r’s are 
real (the absolute value will absorb the imaginary i). Without loss of generality, we may also assume 
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Figure 5. This figure corresponds to the case mi > n 2 . Two consecutive blue edges 
connecting A„ii-i+i to Xj via a /ii, for j<i<n 2 —j + 1, denotes the rewriting of 
the difference Xmi-j+i — Xj as (Ami-j+i + /li) — (Ai + jj,i) The latter is bounded by 
\Xm,i-j+i + Mil + |Ai + Mil and cancelled by the corresponding cross term in (12.301) . We 
perform 712 — 2j + 2 such matching, thus reducing the pre factors in front of Xmi-j+i and 
Xj, in 112.3011 . to ±(mi — n 2 — 1). Similar is the role of the red edges. That is, to write 
Mna-j+i — Mi as (Mna-i+i + Ai) — {fij + Ai), for i = j + 1, ..., n 2 — j + 1 and subsequently 
bound it by |Mn 2 -i+i + Ai| + Imj + Ai| and cancel it with the corresponding cross term in 
(IT^ . Notice that we have performed the matchings in such a way that no edge connecting 
a A and a m is taken twice, i.e. once with a red and once with a blue colour. 


that Ai < A 2 < • • • < Ami that Mi ^ M 2 < ■ • • < Mra 21 that the above expression writes as 


7711 

exp (tt ^(mi - 2j + 1) 
i=i 


n2 _ 

Ami—i+i ^ ^ (*^2 ~ 2j + 1) Mii2—i+1 ~ ^ ^ |Ai + Mil 

i—1 l<i<mi 

i<i <"2 


- 7 r(ni-n 2 ) ^ |Ai| - 71(7712 - wi) ^ IMj 

1<1<7711 


1 )^ 


(2.30) 


Since, in general, mi — 1 > ni — ?72, it is not obvious that the terms (rui — 2j + l)Aj for j = 
1,..., [(n 2 + 1)/2J and j = mi, mi — 1, ...,mi — [(712 + 1)/2J + 1 can be dominated by the terms 
(771 — n 2 )|Aj| and similarly for the corresponding terms involving the m variables, since in general 
ni — 1 > 7712 ~ Till- We will, therefore, need to combine appropriately the Xj and Mi terms in a way 
that it is clear that the combined terms are dominated by the terms |Ai + Mil and thus the factors 
in front of Xj, Mi drop below the values iii — 772 and m2 — mi, respectively (see also FigurejS]). To do 
this we first assume that ttii > 712 . This is no loss of generality, since otherwise we can interchange 
the roles of A’s and m’s in the argument below. However, we still need to distinguish between the 
case mi > 772 and 7711 = 112 • 

In the case mi > 772 we proceed as follows: For j = 1,..., [(112 + 1)/2J we write, by adding and 
subtracting the terms mu for i = j,..., (772 — J + 1) (this is what in Figure [S]we called “matching”). 


{mi - 2j + 1) Ami-i+i - (tui - 2j + 1) Xj 

= (mi - 712 - 1) (Ami-i+i - Xj) + (772 - 2j + 2) Ami-i+i - (712 - 2j + 2) A^ 

»2-i+l ri 2 -j+l 

= {mi —712 — 1) (Ami-i+1 — Aj) + (Ami-i+i + Mi) — (Ai + Mi) 

i=i i=i 

n2-i+l n2-j+l 

< {mi — 712 — 1) (|Ami-i+l| + \Xj\) + E E |A,+M^|. (2.31) 

i=j i^j 
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Similarly, by adding and subtracting the terms Ai, for * = j -|- 1,n 2 — j + 1, we write 

n 2 -j + l n 2 -j+l 

(n2 ~ 2j -|- j+i ~ ~ 2j -|- ^ ^ (m« 2—i-ti ~ ^ ^ 

i=j+l i=j+l 

n 2 -j + l n 2 -j+l 

< ^ \^in2-j+l + ^i\ + ^ lA^j+Ail- (2.32) 

i=j + l i=i+l 

We now write the exponent in (12.3011 as (for convenience we ignore the factors tt) 

mi n2 

^(mi - 2j + 1) + ^(n2 - 2j -f 1) t^n2-j + l E I + Mi I 

j' = l j = l l<2<mi 

-(ni-n2) X! 1-^*1 - (w2-mi) ^ \nj\ 

L(n 2 + 1)/2J mi-L(n 2 + l)/2j 

— ^ ' (m-i 2j -|- 1) (Ami—i-i-i •^i) "f ^ ' (m-i 2j -|- 1) Ami—i+i 

i=i i=L("2+i)/2j-i-i 

L(»2-tl)/2j 

+ ^ ^ (n-2 — 2j -l-1) (/i„2-j+i — /ij) — ^ ^ |Ai + fj,j\ 

j=l l<I<mi 

l<i<"2 

-(ni-n2) ^ \Xi\-{m2-mi) ^ |^j|. (2.33) 

We insert in this the inequalities (I2.31II . (12.321) summed up over j = 1,..., [(n 2 + 1)/2J, and notice 
that since for j = [{n 2 + 1)/2J + 1,..., mi — [{n 2 + 1)/2J it holds that |mi — 2j -|- 1| < mi — n 2 — 1, 
it follows that (12.331) is bounded by 

-(ni-mi-1) ^ |Ai|-(m 2 -mi) ^ 

Since m 2 > mi and ni > mi the desired integrability follows. 

In the case mi = n 2 we do a slightly different matching by writing 

(mi 2j l)(Ami—i+i '^i) (^2 2j -)- l)^/iii2—i+i Mi) 

mi-i rni-j 

= (-^mi—i-l-1 + Mn2—i-l-l) d" ^ ^ ('^mi—i-tl d“ M*) “ ^ ^ ("^i d” Mi) 


n2-j 


i=i+l 

n2-j 


i=i+l 


d- (/r„2-i-i-i d-Ai) -|-Ai) (Mid-Aj) 


i=i + l 


i=i+l 

mi-i 


^ I Ami -i-ti d- /i„2-i-i-i| + E j + 1 H“ I E I Aj Hi I 


n2-i 


i=i+l 
n2-J 


i=i+l 


|Mn2-i-i-i d- Ai| -f I/rj -l- Ai I -l- I/ij -l- 


i=i+l i=i+l 

and finally we proceed as in the previous case. 


□ 


Remark 2.14. The proof of the above theorem and in particular the derivation of relation (12.281) 
shows a way to obtain a contour integral formula for the joint Laplace transform of k point-to-point 
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partition functions via the following recursion (we do not attempt to write the explicit formu¬ 
lae, which would be complicated, or to check the various conditions for the application of the 
Plancherel theorem): 

Let be the partition functions from (1,1) to (mi, ni),..., (wfe, Ufe), respec¬ 

tively, of an (a, Q:)-log-gamma polymer. Assume that mi < m 2 < • • • < m^ and ni > n 2 > ■ ■ ■ > Uk- 
Moreover, without loss of generality, assume that mi < ni. The joint Laplace transform satisfies 
the following recursive relation 






ifc) 


1 

nr(aj + aj) 

i,3 



dA 


.(A) 




r(— 

<mi 


X 


r(Ai-l-Q;j) T{ai + aj)E 

l<z<mi mi + l<z<m2 


e 


— U 2 Z 


(X,c.;d) 
(m 2 , 1 x 2 ) 


UkZ 


(X,c.;a) 1 
(m^, ) 


where = 2,..., k, are the point-to-point partition functions corresponding to a log-gamma 

polymer with parameters (Ai,..., A^^, 0 ^ 1 + 1 , amt; «!, •■■,ani)- 

Remark 2.15. Let us show how our formula for the joint Laplace transform of two equal-time log- 
gamma partition functions provides the corresponding formula for the joint Laplace transform of 
the semi-discrete directed polymer in a Brownian environment studied in |012| and often called the 
O’Connell-Yor polymer. By taking an appropriate scaling limit (this is a straightforward modifica¬ 
tion of Section 4.2 of [COSZ14] in order to include Brownian motions with drift) we can recover the 
semi-discrete directed polymer partition function from the log-gamma polymer partition function. 
More precisely, let be the partition function of a log-gamma polymer at point (mi, Nti) 

with parameters {ai,aj) = (ai,N), for i,j > 1. Let also (i?*(-))i>i be independent Brownian 
motions with drifts ai,i > 1. We then have the following convergence in law, as —>■ 00 , 


log ( 


Nil yiN) 








0<Si<...<Sm,x — 


exp I ^ {B\si)-B\si-i)) jdsi...dsmi_i 


. 2=1 


which is equivalent to 


E 

g “1 

-^ E 




N—¥oo 

. 


where uf = and ti) is the partition function of the O’Connell-Yor semi discrete 

polymer at point (mi,<i). In the case of two points {mi,Nti) and (m 2 , A^t 2 )) such that mi < m 2 
and ti > ^ 2 , we have 


lim E 

N—yoo 



Z 


(N) 

(mi.Wti) 


z 


(JV) 

(ms.JVta) 


= E 




(2.34) 


where and ■ "bhe left hand side corresponds to the joint Laplace transform 

for partition functions at points {mi,Nti) and {m 2 ,Nt 2 ), hence it is given by the formula (12.2111 
with parameters {ai,aj) = {ai,N), i,j > 1. By using the following convergence 


\T{aj+N)) (Af^/v^)-“A 2 


exp(y(Aii-Oj)), 
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Wll -I- W21 
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• • 




Wll -I- Wl2 


Figure 6 . The PNG at time t = 1, 2. The dots represent the entries (wij) of the input 
array, determining the nucleations and the heights of the plateaux. Once created, the 
plateaux expand left/right at unit speed. 


we can take the limit of the left hand side in (I2.34|) and obtain the Laplace transform for two points 
of semi-discrete polymer 


E 




:/ dAs^,(A) n r(-a. + A.) n 0 

/* ^2 —fli , 

._iW2 / 


mi + l<z<m2 
l<z'<m2 


X r(—Ai-h/ii'). 

l<i' <1712 


(2.35) 


3. PNG AND GEOMETRIC PNG 


In this section we define the geometric PNG (gPNG) and use this to obtain alternative integral 
formulae for the joint Laplace transform of the point-to-point partition functions. Naturally, gPNG 
is a geometric lifting of PNG and we start by reviewing the latter. 


3.1. Polynuclear Growth Model and local moves. PNG can be viewed as a con¬ 
struction of an ensemble of nonintersecting paths. In this process plateaux of certain (random) 
heights are created and, once created, they grow horizontally at unit speed. When two growing 
plateaux overlap in the process of growth, the overlapping area drops one level below creating a 
new plateau. At the same time new plateaux may be created on top of the top-most interface and 
the process continues in this fashion. Let us describe this in more detail and let us start with an 
input matrix W = (wij) € (R>o)"^"') which for simplicity we take to be square. This matrix en¬ 
codes the heights of the boxes of unit width that will be created and added on top of the interface. 
It will be clear that PNG can be encoded in terms of local transformations £ij, which we now 
define: For each 2 < i < n and 2 < j < m define a mapping iij which takes as input a matrix 
W = {wij) S (R>o)"^"* and replaces the submatrix 


V ) 

of W by its image under the map 

fa b\ fb Ac — a b \ 

[c d) ^ [ c d + bVc)^ 


(3.1) 


and leaves the other elements unchanged. For 2 < i < n and 2 < j < m, define in to be the 
mapping that replaces the element wn by uii-ip -|- wn and iij to be the mapping that replaces the 
element wij by wij-i + wij. By convention we define £ii to be the identity map. 

Keeping these definitions in mind, we now start describing the geometric construction of PNG 
(we clarify that by “geometric construction” we do not mean the geometric PNG). It is visually more 
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Wsi 


D 


W 22 
' . 


Figure 7. At t = 3, an overlap of height min(/ii 2 ,/ 121 ) — /in drops down and new 
plateaux are created at a; = —2,0, 2. 


natural to consider coordinates (x, t) = {i — j, * + J — 1), 1 < *, J < n, where t stands for time and 
X for the spatial variable. We have: 

• Let us start with a single line of height /i = 0 at time t = 0. At time t = 1 a plateau of 
height ten, of width 1 and centred at zero, is created. It will then start growing leftwards and 
rightwards at a unit speed. 

• At time t = 2 two new plateaux of heights W 12 and W 21 will be created on top of the plateau 
of height wii, which by this moment extends from —3/2 to 3/2. The new plateaux are centred 
at X = — 1 and x = 1, respectively (notice that these are the values of the differences i — j). 
The outer shape of this pile of plateaux will be a line with vertical up and down steps, which 
will have heights 

=wi 2 + wii, for X G [-3/2,-1/2) 

/ij\^ = wii, for X G [—1/2,1/2) 

= W 21 + wii, for X G [1/2,3/2), 


and has zero height otherwise. Subsequently, the new plateaux will start growing at a unit 
speed leftwards and rightwards simultaneously with the underlying plateau of height wn. 

• The third step is a bit more convolved. First, at time t = 3, two new boxes of heights rcia, rcsi 
and unit width will be added on the extended plateau created at time t = 2, above x = — 2 and 
X = 2, respectively, thus creating two new plateaux of heights h^j^2 ^ 21 ^ +'^3i- -^1 the 

same time, the two plateaux which were created at time t = 2, will have grown and they extend 
on top of [—5/2,1/2] and [—1/2, 5/2], therefore creating an overlap, which extends on top of 
[—1/2,1/2]. The overlapping region of height ^ ^ 21 ^ ~ drop below, on the zero 

level, creating a new plateau centred at 0, having width one and height := ^ ^21 ~ • 

At the same time, on top of the topmost plateau which extends above [—1/2,1/2], a new box 
of height W22 will be added, creating eventually a plateau of height = W22 + h^i2 ^ ^21^ • 
So at this stage two lines have been created. The top one with heights 


/i^^^(x, 3) = < 


.(1) 


12 

- 1.(1) 


'•11 

'-21 


— '112 

+ 

^13, 

for 

X 

G 

[-5/2,-3/2) 

= Wii. 

> 

/i^ V /i^V, 

for 

X 

G 

[-3/2,-1/2) 

= W22 

+ 

for 

X 

G 

[-1/2,1/2) 

= Wii. 

> 


for 

X 

G 

[1/2, 3/2) 

II 

+ 

W31, 

for 

X 

G 

[3/2, 5/2) 


and zero elsewhere and the lower one with heights 
/i^^^ (x, 3) = = h \^2 ^ ^'21 ~ 


11 : 


for X G [-1/2,1/2) 


and zero elsewhere. 
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The same expansion, nucleation and drop down process continues. So at time t = 4 the top 
line will have height + rei 3 above interval [—7/2, 5/2). Above interval [—5/2,—3/2) 
a drop down, due to overlap of the expanded plateaux, as well an additional nucleation will 
take place. This action is encoded by the transformation £23 as 


"'12 "13 I ^ ^ 23 ^ / '<■13 /' II ‘22 "12 "13 I 

/122 ^23/ \ h ^22 W23 + h[l^ V h ^22 / 


The lower-right entry encodes the new height of the top line due to the overlap and 
nucleation and the upper-left entry encodes the height of the overlapping area, which will 
then drop down to add on top of the second line . The latter drop down and addition on 
is encoded by the transformation 


(*n . ^*13 ^ ('iit. '^2 + (*13 ^ 


- (.<■')) =: . fig>) 


with , ft. 12 ^) recording the heights of the second line on top of [—1/2,1/2) and [—3/2, —1/2), 

respectively. Similarly, the drop down and nucleation on top of [1/2, 3/2) is described via trans¬ 
formations £32 and £21 as (we show both transformations in one step) 


/ h\V h 


2 

11 

"21 

( 1 ) 

31 


h 


( 2 ) 

12 

( 1 ) 

22 

W 32 


h. 


( 

(1) A l(1) 


A - ^^21 


( 1 ) 


h. 


( 1 ) 

31 


W32 


"12 

^22 

■ V 


•A 




) \ 



) 




/ 

"11 

"12 

...\ 


(i^'S 

•^12 

...\ 



+ {hig ^ - hgi) 

"22 



4 ? 

'^22 





W 32 + hgj V hgi 



h'" 

'^32 



[ 



) 


1 : 


/ 


Finally, a nucleation of a box of height wn will take place on top of the expanded top line. 
This can be encoded via transformation £41 


"31 

W 41 


•^41. 


''-31 

W 41 + 


( ■ 


'-31 

Ml 


V ^ 


The same rules of nucleation and drop down will continue in an inductive manner. This process 
will result in an ensemble of nonintersecting lines. 
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The overall evolution for the first four steps can be given as follows 


fwii Wi 2 Wi 3 Wi 4 ■■■\ 


(h\{ W12 Wi 3 Wi 4 

W21 W22 W23 W24 ■ ■ ■ 

1 .( 1 ) 

W21 W22 W23 W24 ■ ■ ■ 

W31 W32 W33 W34 ■ ■ ■ 

W31 W32 W33 W34 ■ ■ ■ 

W41 W42 W43 W44 ■ ■ ■ 


W41 W42 W43 W44 ■ ■ ■ 

\\ \ \ \ ) 


\ \ y 


"'ll "12 


W 12 , Wi4 

f,(l) 

"21 '^22 W23 W24 

W3I W32 W33 W34 

U>41 U ;42 UI43 W44 


//.S? 

( 2 ) 


..A 




^41 O (^210-^32) O (^120'^23)0'^14 


^21 

^31 

Ml 


'M2 

'M2 

'M2 

^42 


■^ 310 ^ 220 ^ 13 , 




'n3 'n4 

uW 

"23 '“^24 

W 33 W 34 

W43 W44 


(1) ;,(!) 


h^^'> 

"21 

h^^'> 

"31 


,(i) 


"-22 

W32 W33 


13 '*^14 

W23 W2i 




WAl W32 W33 


W34 

W44 


...\ 


Introducing the transformations 


I ^l,j — i-\-l 0***0 — \ O ^ ^ j 

I j+ 1,1 o • • • o £ 2 —i,j —1 o £ij , i ^ j-i 


(3.2) 

we see that, given a nucleation matrix W = (wij) G (lR>o)"^", we can write the PNG evolution in 


terms of transformations r® as 


< o «-i °r^ {ri o---orl^) o-- - o (rlorl) o rl (W). 


(3.3) 


The output matrix 'H will encode the heights of the line ensemble. Let us, also, remark that the 
nucleation matrix W = {wij) need not actually be a matrix, but it can rather have any shape. An 
example with which we will be occupied is when W = {wij , l<i+j — l<n), corresponding to a 
triangular array. 


3.2. Geometric png and local moves. Rewriting the transformations £ij in definition 
of p.l|) by replacing (+,V) with (x,+), i.e. performing a geometric lifting of the formulae therein, 
we clearly arrive to the geometric local transformations lij defined in (EZD. Similarly, writing the 
transformations r* in (13. 2p in geometric form we arrive to the transformations of (12.101) . The 
geometric version of matrix equation (ESI, i.e. 

n = gPNG(lT) := £>” O o Q^-^) o • • ■ o {qI o ■ ■ ■ o q\) o • • ■ o {gl o gj) o p) (IT) 

can be taken to be the definition of the Geometric PNG corresponding to matrices. Notice that even 
though gPNG is a composition of the same local moves gj and kj employed by gRSK, the sequence 
of the compositions is different. However, it turns out that gPNG(lT) = gRSK(lT), whenever W is 
a matrix array. 

A reason why we are interested in the gPNG formulation is given by the following proposition, 
which says that when gPNG is applied to a triangular array W = {wij : i,j > 1, l<i+j — l<n), 
the elements hij, with i + j = n + 1, of the output array are the point-to-point polymer partition 
function from (1,1) to (*, j), i + j = n -|- 1. In the case of a triangular array IT = {wij '■ i,j > 1, 1 < 
i -|-j — 1 < n) the output array via gPNG is a triangular array 'H = {hij : + j — 

defined by the sequence of local moves as 

n = {glo---o g^)o---o{glo gl)o g\{W) 


(3.4) 
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The properties of gPNG on triangular arrays as defined in (13.41) will be the object of the following 
propositions and will eventually lead to an alternative formula for the Laplace transform of the 
joint distribution of the point-to-point partition functions at a fixed time (cf. Theorem l3.5|l . 

Proposition 3.1. Let W = Wn '■= {wij '■ i,j > 1, 1 < * + J — 1 < n) be a triangular array and 
'H = {hij : *, j > 1, 1 < j j — 1 < n) = gPNG{W) = o • • • o g”) o ■ ■ ■ o o g^) o gl{W) be the 
geometric PNG array. Then 

hij = Zij = Wij, for all {i,j) with i + j = n + 1, 

TTGlIij (ij)G7r 

where Ily is the set of directed paths (in matrix notation directed means down-right) starting at 
(1,1) and ending at (i,j). 


Proof. We can recursively write gPNG in 

the form 







( 








■Hn-1 

W2,n-1 

Tin ■■= gPNG(W„) 

II 

o 

0 q\ 

\Wnl 

Wn-1,2 




( h'li 

h'i2 



Wln\ 



h '21 

h' 

'^22 


W2,n-1 



II 

o 

o 








h' 

'to- 1,1 

wn- 1,2 





\ Wnl j 


Wln\ 


/ 


(3.5) 


Notice that the transformations g) only change the entries on the diagonals (i, j), (i — 1, j — 1),.... 
Moreover, if hjj = Zk^n-k, for k = 1,2..., then transformation g^_f.j.i will produce the entry 
hk,n-k-\-i = Wk,n-k+i{h'k_i^n-k+i + which by the basic recursion of polymer partition 

functions equals Zk^n-k+i- D 


Proposition 3.2. The gPNG, in logarithm variables, is volume preserving, that is, ifTl = {hij, 1 < 
i-\-j—l<n) = gPNG{{wij)i<ij-j-i<n), then the map 

{\ogWij, l<i+j-l<n)H> (log hij, l<i+j-l<n), 

has Jacobian equal to ±1. 


Proof. It follows immediately from the representation of gPNG in terms of transformations g) and 
their volume preserving properties. □ 

In the setting of gRSK the energy of a matrix X G (lR>o)"^’”, was defined as 


£{X) 



Xij 


(3.6) 


where the summation is over I<i<?T,, 1 < j < m with the conventions Xij = 0 for i = 0 
OT j = 0. We can extend this definition for general arrays, with the convention that Xij = 0 if 
{i,j) does not belong in the index set ind{X) of this array and in particular for triangular arrays 
X = {xij, l<i+j — l<u). Similarly to part 2. of Theorem 12. II we have 


Proposition 3.3. Let W 

Then 


{wij : 1 < * + j — 1 < n) be a triangular array and TL = gPNG{W). 


£{n) = ^ 


1 
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Proof. The proof is similar to the one of Theorem 3.2 of |OSZ14| and it is based on the invariance 
of the energy under transformations p*. For a triangular array X = (xij , * + J < n + 1) we define 


= — + E 




Xii 


E 


i-j-j—n-j-l, i<k 
or 


i-i-j—n-j-l, i>k 


Clearly, = £(X) and £^(X) = ^Iwij. Let us assume that the 

claim is valid for n — 1, i.e. £{{xij)i<i+j<n) = X]i<i+j<n Using the recursive structure of 

gPNG, cf (13.51) . it is enough to show that 


£t^{gt.k+iiX)) = £!:(X), for k = 


(3.7) 


Let us denote by X' = (x£)i<i<j<n = Pn-k+iiX)- When fc = 1 we have that x£ = Xij for all 
ihj) ^ (1)''^') and x'l^ = winXi^n-i- Therefore, 


^hix) — -h ^ 




Xll 

1 




i+j—n+l,i>l a 


1 

Wii 


Xi-lJ + XiJ-l ^ 1 ^ 1 

a;ii Xii Win .. . 


= —+ E 

nr-t 1 

z+j<n 


i>2 


1 


= —+ E 

7*1 1 


H” — l ^l,n—1 


1 


E — 


3/11 ^ 1/n 


= £ni 9 UX)). 

Let us now check (liTTll for 2 < k and k ^ L^/2J +1, if n is odd. It suffices to show that for 


y^'/ 1,3-1 ^'i-i,3-i \ 

V X£ X', , 1 X', n / 

•'J 

E"( 


ij-1 






^fc,n—fe+1 


^i,j — 1 


— 1 ,j 


where summation is over (i, j) = (fc — r, n — fc + 1 — r), with 0<r<fcA(n — fc + l) + l and 
is similar except that l<r<A:A(n — fc+l) + l. This is confirmed by the relations, valid for 
ihj) = (k — r, n — k + 1 — r), r = l,...,fcA(n — fc + l) + l 


—1,3 ^ Xi^j — 1 X'l —1,3 I Xz^j — 1 


1 


1 


*,3 


*.3 




+ 


1 




1 




X' 


k,n—k-\-l 


^k,n—k-\-l 


which follow from (12.111) and the fact that x^i j = a^iii.j and x' = Xij±i. In the case that n is 
odd and k = [n/2j +1, we use the same set of relations and in addition that 


‘-11 




1/1 


+ 




2^11 V2;12 2:21/ 


□ 


The following proposition is the analogue of the fact that geometric RSK preserves the (geometric) 
type of an input matrix cf. ProDOsitioi l2.7l 
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Proposition 3.4. Let W = be a nucleation array and 'H = {hij)i<i+j-i<n = 

gPNG(W). Then for every p, q such that p + q = n + 1 or p + q = n it holds that 

pAg—1 

hp—r^q—r — '^ij ; (^■^) 

r—0 

where we notice that the rightmost product is over all the input elements Wij inside the rectangle 
with downright corner (p, q). 


Proof. The proof goes by induction, so let us assume it is valid for n — 1. Let 

T~L ■= i<n—1 = gPNG((rcjj )i<2_|_j_i<yj_i) 


then 


'H = Qfog'f 1 o . . . o o g^ 


= 01° Q2 ^ 


' o qI-i o 0 }, 


/ 

n' 


wn- 1,2 

\Wnl 


/i'll 

h'i2 

h '21 

h' 

''-22 

h' 

''•n-1,1 

Wn-1,2 


yjln\ 


W2,n-1 


/ 

W2,n-1 




\ Wnl 

Notice that the transformations only change the entries on the diagonals 

{(bj)' J ~ * = n — 2k + 1} with k = while leaving invariant the in between diagonals 

{(i,j):j — i = n — 2k} with fc = 1,n — 1. For this reason we have that hp-r,q-r = ^p_r q-r 
p, q such that p + q = n and r > 0. Therefore, (13.81) is immediately satisfied since by the inductive 
hypothesis 

pAq—1 

= n for P + q = n. 

r—0 

Let us now check dim when p, q are such that p + q = n + 1. We consider the product 

pAq—1 pAq—1 

hp—r,q—r — ^p,q r,g—r 

r=0 r—1 

Using the definition of g}, cf. (12.101) and (12.111) . we can write the above as 


pAg —1 

'^p,g(^p-l,g + ^p,g-l) 


1 


h' h' 

"'p—r+l,q—r "'p—r,q— r+l , L' I 

\'^p—r — l.q—r ' ^p—r,q—r—l)’ 


-- h' h' h' ^ ^ 

— ^ p—r^q—r p—r-\-l,q—r ' p—r,q—r-\-l 

The terms j and h^_i j + will telescope and so the above equals 


Wr, 


pAg-l 7 / rrpAg-l 7 / 

r—1 ^p—r-\-l^q—r llr=l '^p—r,g—r +1 


nPAg-i 7 / 
llr=l p—r,q—r 

Using the inductive hypothesis this equals 

ni<p,i<g-l ni<p-lj<g 




n 


i<p-l.i<g-l 


n 




and this proves the claim. 


□ 
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Combining the above propositions we arrive easily at the following theorem, which gives an 
alternative formula for the joint law of all the point-to-point partition functions at a fixed time n 
of a directed polymer with log-gamma disorder: 


Theorem 3.5. We consider the log-gamma measure on arrays W = (wij, 1 < — l < n) defined 

by 

P(dtT):= Yl r(a,+ d ,7 > 0, (3.9) 


The joint distribution of the log-gamma polymer partition functions at time n 

Zji .— ^Zp qj p -\- g — n -\- , 

is given by 


F{Zn G da;) 


X 



1 

n*,j r(ai -f dj) 


n 


dxj 

Xi 


( Ylj-i=n-2k+l 
Y[j-i=n-2k 




n 


1 

rijj r(ai -I- dij) 


i=l 


dxj 

Xi 


f T{j-i=n-2k+l \ JJ- dUj 

\nj-i=n-2fe+2 / i+j<n 


(3.10) 


where T(a;) := [{tij), 1 < i < j < n with tij = Xi when i -\- j = n -|- l) and we also follow the 
convention in the above formula that if {i,j) do not satisfy l<i-\-j — l<n, then tij = 1 . 


Proof. Proposition 13.31 shows that J^i j — ^(T), where T = T(a;) when tij are set equal to Xi 

for i-\-j = n -f 1. Proposition 13.21 shows that ])[^ ^ dwij/Wij = Yiij dUj jtij. Finally, from Proposition 
1 shows that 


n 




n 


j — i—n—2k+l 


n 


j—i—n—2k 


and 


n—k-\-l 

n 

i=i 


Wk,i = 


n j—i—n— 


j-i=n-2k+l 


n j—i—n— 


2k+2 


which raising to the powers —d„_/c+i and —otk, respectively, gives the corresponding factors in (I3.10|l . 
Putting the pieces together provides the push forward of the log-gamma measure on the gPNG array, 
which when integrated, fixing the entries tij = Xi for i-\- j = n -I- 1 leads to distribution (13.101) . The 
proof is now completed by recalling Proposition 13.II which identifies the entries {ty : i j = n1} 
of the gPNG array with the point-to-point partition functions of the directed polymer. □ 


Remark 3.6. It is worth noting that the function (xi,..., a:^) appears to have an integral 
structure of the same flavour as that of Whittaker functions, however the structure of the potential 
f(T(a;)) is different. This is indicated in Figure |S] where the convolutions in E{J{x)) are depicted 
(compare also with the analogous diagram for Whittaker functions in Figure [TJ) 


4. Towards the Airy process 

In this section we formally check the validity of the following conjecture 

Conjecture 4.1. Let 

{mi,ni) = {N -tiN^/^,N + tiN^/^) and (m 2 , n 2 ) = (A-f A - (4.1) 

with ti,t 2 > 0, and ^(m 2 ,n 2 ) point-to-point partition functions of a —log-gamma 

polymer (i.e. an {a, d)-log-gamma polymer with ai = 0 and dj = 7 for i,j >1/ If Ai{-) denotes 
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Figure 8. This figure shows the structure of the potential £{j{x)) appearing in the 
function (ii,in) in (13.1011 . The arrows in the picture show that the bonded vari¬ 
ables are convoluted in the integrals, i.e. there are exponential weights exp(— 
or exp(—tij/tij+i). The pointer of each arrow indicates the variable that lies in the de¬ 
nominator in the previous fractions. The arrow pointing towards tn corresponds to the 
term exp(— l/tn). 


the Airy process (which will be formally defined in Section o and the functions -F, G are defined 
as 


G'(z) ;= logr(z) - logr (7 - z) +/-yZ and F(z) ;= logr(z) + logr (7 - z), 
with = — 2 T'( 7 / 2 ) and 111 ( 2 ;) = [logr]'( 2 ;) the digamma function, then 


(4.2) 


log ^(mi ,ni) - Nfj log -Nf^ 

(c7)-i1V1/3 ’ (c7)-17V1/3 

where o' - i I-l/3 7 _ c7(F"(7/2)) 

wnere Ci — { 2 ) ? ^2 ~ 2 G'"( 7 / 2 ) 


(d) 


N—¥oo 


> {Ai{-cjti) - cl tf, Ai{clt 2 ) - cl 4) 


and ( 


F"h/2) 

■G'"(7/2) 


Even if formal, our check of the validity of the conjecture makes some non trivial steps. In 
particular, we re-express formula (12.201) as finite series related to the Fredholm expansions of the 
Laplace transform of a single partition function and then we reveal a structure of block determinants 
that relate directly to the Fredholm expansion of the two-point distribution for the Airy process. We 
can then pass to the limit on each individual term in this expansion and establish that it converges 
to the corresponding term in the Fredholm expansion for the Airy process, cf. Proposition 14.101 
However, what makes our computation formal is that we were not able to justify that this limiting 
procedure can be done uniformly for all terms in the expansion. The difficulty arises from the 
presence of a cross term 114.2211 involving gamma functions. A naive estimate of this cross term with 
a uniform, constant bound on its modulus, when the values of the arguments are in a compact set, 
produces a bound of the form C"™, which cannot be compensated by the factorial l/n!m! in the 
expansion to truncate the series. Besides such an extension of Proposition 14.101 to hold uniformly 
to all orders in the series and (for the same reason) a formal interchange of integration in relation 
(I4.12|l (which, however, can be made rigorous in the case of the O’Connell-Yor polymer), the rest 
of our computations are rigorous. 


4.1. Auxiliary results. The Airy process Ai{-) is a stationary, continuous process with finite 
dimensional distributions given by 

P(Ai(ti) < fi,...,Ai{tk) < Cfc) =det(/-fAif)i2(p^^..._t,,}xiR), (4.3) 

where for ^ = 1 ,..., k we define f(t^, a;) := and Ai is the extended Airy kernel defined by 

f g-Ht-t )Ai(^ + X)Ai(^' + A)dA, if t > t', 

Ai(t,e;t',e') := < 

[ - /°oo Az(e + \)Ai{C + A)dA, 


A t < t', 
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In the above integral Ai{-) is the Airy function, which has the integral representation 

Ai{x) = ^ ! e^-^^dz (4.4) 

3 

The contour of integration consists of two straight rays, one starting at infinity and ending at the 
origin in an angle —tt/S and the second departing from the origin at an angle tt/S and going to 
infinity. Let us mention that the Airy process can be viewed as the top line of an ensemble of non 
intersecting paths, which is known as the Airy line ensemble, whose statistics are also governed by 
kernel Ai. 

The determinant in (14.31) is a Fredholm determinant on ..., tk} x R). Let us recall that for 

a trace-class operator K acting on L^(A’, p) as Kf{x) = K{x,y)f{y) d^{y), for some measure 
space the Fredholm determinant is defined as the series 

det(/-f A:)l 2 (a') := 1 + X! “1 / det (A:(a;i, dp(a;i) • • • dp(a;„). 

In particular, the Fredholm determinant in (14.31) can be written as 


det(/ - f Ai f)L 2 ({ii,...,tfc}xR) 
" (- 1 ) 
n—1 




E 


/ OO nOC 

■■■ det((fAif)(s„2/*;Sj,%))„^^d2/i---dy„ (4.5) 

-OO J —OO 


It will be useful to rewrite the above formula in terms of a block determinant. We do this as 
follows: We divide the variables si,..., s„ into k groups according to which value out of {ti,..., tk} 
they take. The number of ways we can do this assignment so that for j = 1,..., fc we have Uj variables 
taking the value tj equals 

/ n \ n\ 

Taking this into account and the fact that n = ni -f • • • -I- rife, we can rewrite (113 as 


^ -hnfc 

det(/-fAif)^2(p^_ = 1+ E 

lii. Ilk- 

ni ,...,rtfc = l 

■ OO / \ "iH l-"fc 

det Y[ dyj, 

\ + j=i 

where n,. = the determinant that appears above is a block determinant. 

Let us record the following Lemma, which is a generalisation of Lemma 4.1.39 in |BC14) . The 
latter has formalised an idea, which has appeared in the physics literature earlier, e.g. [CLeDRiO] . 
and reduces the study of the distributional limit of random variables to the asymptotics of the 
corresponding Laplace transform. The proof of Lemma 14.21 is a straightforward generalisation of 
that in |BC14| and so we omit the proof. 



Lemma 4.2. Consider a sequence of functions {fn}n>i mapping IR —>■ [0,1] such that for each n, 
fnix) is strictly decreasing in x, tends to one as x tends to —oo and zero as x tends to oo and for 
each (5 > 0 converges uniformly to l{a;<o} Define also ff{x) = fn{x — r). Consider two 

sequences of random variable (A„) and (y„) such that for each ri,r 2 G R, 

E[/;n^n)/;ni"«)]^p(ri,r2), 

and assume that p{ri,r 2 ) is a continuous joint distribution function of two random variables. Then 
(Xn,Yn) converges weakly to a random vector {X,Y) which is distributed according to P(A < 
ri,Y < r 2 ) = p{ri,r 2 ). 
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We will apply this lemma by considering the functions fN^x) = e ' 
criteria, and choosing the parameters of the joint Laplace transform as 

and 112 = 




which satisfies the 


111 = 


(4.7) 


where = —2'I'(7/2) and 41 ( 0 ) = [logr]'(z) is the digamma function. For i = 1, 2, we have 


— UiZf. 




N 


WT^ 


By Lemma 03 for each pair (ri,r 2 ) G it suffices to compute the limit 


lim E 

N—^oo 


IN 


ivV3 


IN 


log ^(m2 ,n2) -Nf^ 
iVl/3 


= P-i{ri,r 2 ) 


and show that this coincides with the joint distribution of the Airy process at the corresponding 
points. In other words we will need to compute the limit of the Laplace transform E 
with the parameters set as above. 

We will make use of the following theorem due to Borodin-Corwin-Remenik |BCR.13| . We will 
also use the notation: iL£s is a vertical line in the complex plane with real part equal to ±5, 
is a horizontal line segment from — -I- lM to 82 + M and is a circle centred at zero with radius 

S. 


Theorem 4.3 ( [BCR18) L Fix 0 < (52 < 1,0 < (5i < min{(52, 1 — 152} and oi, .., oat such that \ai\ < (5i 
and let F be a meromorphic function, which is non-zero along and inside Cs^ and all of its poles 
have real part strictly larger than 82 ■ Moreover, assume that for all n > 0 


J±e,„ Jt' . \M\^oo 


^ ±^ 5 , 

Then 


f ^ ^ F( ■) 

/ dwsN{w) T{af - Wj)Y[-^;p-=det{I-G K)l 2 (^Cs,)- 


The kernel K : L‘^{Csf) —^ L‘^{Cs,f) in the above Fredholm determinant is given by 


K(v, v') 



dw TT F{w) r(i; 

w — v' sin(7r(i; — w)) F{v) r(ii; 


Of)' 


In particular, recalling Theorem \2.1(k we have that the Laplace transform of an {a, d)—log-gamma 
polymer partition function is expressed as 


E 



det(/ -I- Ku)l^(Cs^), 


(4.8) 


where the kernel : L‘^(Cs,^) —>■ Lf{Csf) equals 


Ku{v,v') 


1 f dw TT u™I-ff{—w) -A- r(i;-|-Q:j) 

2ttl Ji,. w — v'sm{Tr{v — w)) -IJ-r(ii;-|-d^-) ’ 

2 J — i 


where notation F(),(-) was introduced in DeRnition \2.!A 


Remark 4.4. Even though, in general, a Fredholm determinant is given by an infinite series, an 
inspection of the proof of the above theorem [BCR13] shows that (14. 8|) is finite rank and the terms 
beyond order N in the series are zero. 


2>'^2) 
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4.2. Road map. Let us describe the steps we make towards Claim 1. First, in Proposition (14.51) 
we use Theorem l4.3l in an iterative manner, in order to rewrite formula (|2.20l) as a double, finite series 
of terms related to expansions of Fredholm determinants like (14.81) (denote by Im)i the generic term 
in these double series ). We note that the proof of this proposition is formal, as it contains a formal 
interchange of integrals in (|4.12l) . On the other hand we also note that this formal interchange of 
integrals can become rigorous, in the case of the O’Connell-Yor polymer, cf. Remark 14.61 However, 
for the sake of exposition and uniformity of the text, we have preferred to present the calculations 
in the case of the log-gamma polymer, even if these contain a formal step. 

Then in Lemma we make an observation that certain block determinants arise in the above 
mentioned expansion, which point towards the expansion dMl) of the Airy process. 

In Proposition 14.81 we further rewrite our formula, so that it is in a form more suitable for 
asymptotics. The asymptotics are then carried out with the use of steepest descent method. For 
this we need to make suitable contour deformations and also estimate our integrands along the 
deformed contours, cf Lemma 14.91 and (14.171) . 

The last step we perform is to show that using the steepest descent method, any fixed term in 
our double expansion, denoted by Im,rn converges to 


(- 1 ) 


n+m 


n! TO!(27r6)’”+"' 


pOQ p' 

</c^r 2 +C 2 i 2 


clr 2 +c^t^ 


i=i 


n+m 


n 


c/j-i-l-cRf Jcln+citf j^n+l 


(4.9) 


/ 

det 

V 







Ui+n j 



/ 


and when this is summed over all possible n,m> 0 it coincides with the probability 


¥(^Ai{-clti) < cjri + c^tl ,Ai{cJt2) < c +2 + cjtl ^ . 

Then the use of Lemma |4.2I would confirm Claim 1, assuming one could obtain a good control on 
the tails of the series involving the terms 1 ^ 2 . 


4.3. Main computations. We now follow the roadmap and start by rewriting the contour 
integral formula provided by Theorem 12.121 for the two point function. 


Proposition 4.5. [modulo justification of formula (14.121) ] Let Z(^rni,ni)y ^(^^ 2 , 112 ) partition 

functions from (1,1) to (mi,ni) and (rn 2 ,n 2 ), respectively, of a (0,"f)-log-gamma polymer. Fix 
0 < ^ < 7/2 and 0 < < min{^, 1 — ^}. We assume that mi < ni, m 2 > n 2 , m 2 > mi and n 2 < ni. 

Then the joint Laplaee transform is written in the form 


E 




n 2 7n ~\ C C 1 1 

= ;—u-—^— 1 — / dvdw / dvdw det ( - 2 det ( - 2 

^ n!m!(27r(.)"*+" 7 J \Wk-ViJnxn Vw^-v^/ 

n—0 m=0 

sin(7r(ufe - Wfe)) M2'' \r{wk) J 

^ A ^ M^ / r(7-w,) yY^A)A^ 

sin( 7 r(v^ - w^)) u\‘ V r (7 - v^) J ^(w^) J 


X 


A A r(7 -We- Wk)T{j - Vg 
r(7 - wf - Mfe) r(7 - v<> 


Wk) 


where the first set of contour integrals above are over (Cs^ x £ 5 )” 


and the second over (Cs^ x £ 5 )”*. 
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Remark 4.6. The two-point function of the O’Connell-Yor polymer in Remark 12.151 can also be 
written in terms a Fredholm-like expansion: 


E 


^ — UiZ — 


-EE 

n—0 m—0 


1 


n!m!( 27 rr)"‘+" 


/ dvdw [ dvdw det ( -- '] det ( -) 

J J nxn — \/£ / mxm 


n 


sin( 7 r(ufc - Wk)) (ui/u 2 )’''“ V J 


Wk — Vi / nxn 

m2 —mi 

exp I 


k— 


n 


exp(^K-u,)) 


sin( 7 r(vr - w^)) \r{wt) J 

^ TT TT r(?ilfc — Vill)T(vk — Ml) 

M 1=1 


(nMl)Y^ (t2, 2 20 

lrta;T 


where the first set of contour integrals above are over {Cs^ x ig ')” and the second over x tg/ )"* 
with (52 < (5^ < (52 < i^i- The proof follows the same steps as in the upcoming for log gamma 
polymer. In this case, however, the interchange of integrals can be easily justified by Fubini, due to 
the super-exponential decay of the terms exp and exp ~ )) along vertical 


axes. 


Proof of Proposition 14.51 . Let us start with a general (o;, d:)—log-gamma polymer such that 
jofil < (5i and \aj — ^\ < (5i and with formula (I2.2()F We multiply and divide inside the integrals of 
(I 2 . 20 |l by Y\T=i nj=i = 0^=1 r("^i + %) and rearrange the terms to write ( 12.2011 as 


u7^' F° 


„ MI TT T^C „ I \ 1 TT “1 ' n2,ni{^i) T-j- 11^ = 1 + %) 

j dA Srn,i (A) II r( 0 ;i/H-Ai)|| -ai r-n { \ 11 TT ^2 p/ i ^ ^ 


/ dp s„jp) Jl r(-dj/-f pj)P[ 

(■^5 + t) ^ ^^jij''^'>T'2 J —1 




mi ,m2 


(dj) t 4 nil\ + dj) 


n 


Fmi.m2(%) j=l ni=l 


dA s„ 






((A) n r(-a.,+A.)n E. 


ni V 
k-j pA 


1 / \ TT t^/ ~ I I TT “2 ' mi ,7712(^1) 

I dps„Ap) 11 r(-aT^a--—— - —■ 

(^5+-,)^2 i —1 ^2 ^mi )^mi ,m2 ) 


(4.10) 


We will now rewrite the integral over p using Theorem 14.31 in terms of a Fredholm determi¬ 
nant. Observe that this integral represents the Laplace transform of the partition function of a 
(Ai,..., Xmi, Omi+i, Q:m 21 di,..., d^^)-log-gamma polymer on [1, m 2 ] x [1, 77 , 2 ]. Therefore, the inte¬ 
gral will be the same if the parameters are changed to (Ai -|- 7 ,..., Am^ -I- 7 , Omi+i + 7 , 0 'm 2 + 

7 ; di — 7 ,..., d „2 — 7 ). This will be technically convenient later on. We then have that 


-f^i vzX 


'(d+7)"2 


dp s„ 


(d) n +H)Y[ 


^mi iP'j ) ^mi ,m2 (dj ) 
L(ai)Fmi.m 2 («j) 




^2^^'Fmr(M7)F))+7m2(A^7) 

=' 1 “^“'’'^^ FAAaj)F“i,„i2(“7) 


/ dps„2(p) P[ r(-dj/-f7 

l<j,j'<n 2 

[ H 1 1 TT FI" " 4- N t4 “ 2 ^ FA7"’^(~di)FA7)m2(~di) 
/ dp s„Ap) M r(-ay-f7-pA M — 

!<.,/<». i= "■ 

^A 


i=1^2“"’^^ F^A'A7)FAi.m2(«j) 


— det(/ -I- 
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where the kernel Ku^ '■ L^{Csi) equals 

dw 7T 




w-v' Sin(7r(u - w)) u1^ f^+y r(u; + dj - 7 ) ‘ 


n 


r(u + dj — 7 ) 


Notice that all the poles in the product F^+'''(—(—have real part strictly larger than 6. 
This is a consequence of the assumption (5i < ^ < 7/2 and \ai\ < 5i. Moreover, the decay conditions 
in Theorem 03] is verihed by using the asymptotics of gamma function 12.151 Inserting this Fredholm 
determinant into (|4.10|) we obtain that the latter equals 
„ mi _ -Ai 


/ dA (A) n r(-a,+A.) n 


det(/ + K^Jl2^Csi) 


(4.11) 


The next step is to expand the Fredholm determinant in (14.1111 and then expand the determinants 
involved in this expansion. Recalling Remark 14.41 which implies that the series expansion of the 
Fredholm determinant contain only 712 number of terms, we have 

n 2 


' 1 f 

det(/ + K^^)mcs,) = det{K^^{vk,ve))nxn 

71 

dv n 


n—0 

V V 

n' 

ra=0(TeS„ 

Al2 




/c=l 


= EE 

n—0 (T^Sn 


-E 


'i!(27rt)’^ 


dv / dw J|- . . E -^ 

(Cjp)" Wk - fa(fc) sin(7r(wfe - Wk)) 

{-Vk)^^kVm2{-Vk) 


Vk 11^+7 
^2 ' '^1 


n—0 


d(27rt)’^ 


dv 


dw det 


'(Cdi)" 


'(dY 


1 " 

(^—) n 

\Wk — Vl/ nxn 


ve/nxn sin(7r(-(;fe - Wk)) 


where in the last step we used the fact that 


CT^Sj 


(-irn 


1 


(-r;fe)F“+;i,(-Ufe) Ft-'^iwkY 

(; 


,Vk p '*'+7 

^2 ' 


= det 


1 


V Wk — V£y nxn 


k^l '^k Vcr(k) 

Inserting the Fredholm expansion into (14.1111 we write the latter as 
1 r . -r-r . . -rT u 

_ . \ rj / 


n 2 


n!(27rt)” ' 


n—0 


(A) 


TT , \ A TT “1 F“ (Ai) 

,2'<mi i—1 1 niv 


/ dv / dw det 

liCs.Y hdY 


1 


n 

\Wk — Vp y nxn 


^^kV{-Wk)FYY?m2^-Wk) K:^{vk) 
Wk-Vt)nxu ^^^sin(TT{vk-Wk)) uY {-Vk)FYit?m 2 {-Vk) KY^iwk) 

rearranging formally the integrals (however this rearrangement can be justified when 2 n < ni — mi) 
and using that ]Xk=i ^mVi~'^k) = OLi Y\T=i F(Aa + 7 - ^fe) = OE F;(“'"(Ai) and similarly for 
nLi Fi);tn —Vk) we write the above as 

_ _F“+Xi,(-Wfe) F“-T'(z;fe) 

. . f I \/ f f I \A/ f I l-> ■ ■ I ■ ■ 

l"o^!( 2 ^ 0 ”Edi)' 


n2 

E 


dv 


dw det 


'(dY 


(—^) n 

\Wk — Vpj nxn 


I'fc pQ^+7 


mi ,m2 


I dA Suii (A) 


xn sin(7r(z;fc - Wfc)) u. 

;,i'<mi i=l ^1 (A, 


{-Vk) K 2 "^{Wk) 


F“(Ai)Fr“(Ai) 


(4.12) 
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For A G C, set 

. F“(A)Fr“'(A) 

•“n 

Multiply and divide the last integrand by nlL\ F“+'^(--yfe)/F“+T'(-'u;fc) 

and rewrite the above as 


E 


dv 


dw det 


1 /<- 

(^—) n 

V Wh — Vf / nxn . 


F“+^(-R;fe) FtTHvk) 


^ n!(27rt)" 7(c,,)" V)" '^Wk - veJnxn sin(7r(pfe - u>fc)) it“'= F^p{-Vk) F“, '^{wk) 

-Xi r /AX 

dAs™,(A) n r(-a.>+A.)n ■ ( 4 . 13 ) 






We will now use Theorem 14.31 to compute the integral with respect to the A variables as a 
Fredholm determinant, i.e. 

p 777-1 —Ai £ / \ ^ 

/ dAs^,(A) n r(-a.. + A.) n ""ia. ^ = det (J + 

his)-- 

where the kernel ^ L^{Cs^) is given by 


T^n,w,v 

Ku’ 


1 


(v,v') = :^ 


dw 


i™C.u-,«(-w) ^ r(v-l-aj) 


27rt w — v' sin(7r(v — w)) r(w -|- Oj) 


n 


(4.14) 


All the conditions of Theorem l4.3l are verified: The poles of the function fn,!«,?;(—A) are dj and ■j — w 
which have real part strictly larger than S. The decay conditions are verified by usual estimation of 
gamma functions. Formula (14.131) is then written as 


E 


1 


. / dv / dw det 

^ n!(27ri)" 7(d)- \Wk-ve/nxn 


n 


(^) 

\Wk- VI J1 

uT K-'<{vk) 


sin( 7 r(r;fe - Wk)) ul’‘ FZV{-Vk) F“, '^{wk) 


det + 


L^(Cs,) 


We now specialise and restrict attention to a (0, 7 )-log-gamma polymer, i.e. when Oi = 0 for 
i = 1,2,... and dj = 7 > 0 for j = 1,2,.... Moreover, we expand the Fredholm determinant 
det (/ + J and perform the same calculation as above to arrive to the result. □ 

The next lemma is crucial as it reveals the block determinant structure, which relates to the 
Fredholm expansion of the two-point distribution of the Airy process. 

Lemma 4.7. For any sequence of complex numbers {wi)i>i, (r’i)i>i, (wi)i>i, (vi)i>i, such that for 
some J > 0 it holds that 5R(wi), 5R(wi) < 7/2 and Ji(rii), Ji(vi) < 7/2 fori = 1 , 2 ,... and furthermore 
— Vj) >0 and 5R(wi — Vj) > 0, the following identity is true 


det 


(7 - Wj; - ZIfc ) (7 - V/ - Wk) 

Wk - VeJ nxn'^'' \wk - yeJ mxm {j - - Wk)i"t - - Vk) 


^ - n m 

(^—) det (^—) n n 

\Wk- — Vpynxn VWl — V/? / 777 xtti 


(4.15) 






f (0,oo)” 


dTi • • • dXn+m det I 


)) (-( 

/ 77 X 77 V 


-Xj + „(j-Wi-Wj 






4, 


A 


Proof. Expanding the Cauchy determinants in the left hand side as 

ni<fe</<„(w'fc - wi)('>ii - vk) 


det 


(^) 

\Wh — VP y 


Wk — Vi y nxn 


Ul<kJ<n(.Wk - Ve) 
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and recombining with the product term therein we can write the left hand side of (I4.15|l as a block 
Cauchy determinant 


det 


(^) ( ^ 1 

V Wi—Vi ) V lUi+Wi—7 / 

\ ^/nxn \ ' / 7 

i^) (^) 

\ ■y—Vi — v-j / \ Wi—v-j / 

^ / mxn \* ^ / m' 


= det 




nxn 

(r (/“ e-(--''^)dx) 


ddx) 


nxm 
mxm 


Notice that the condition 5i(wi) < 7/2 and Ji(vi) < 7/2 for i = 1,2,..., ensures that 

the above integrals are well defined. The final step is to use the multilinearity of the determinant, 
in order to pull the integrals out of the block determinant and hence arrive to the right hand side 
of giSD. □ 

Using Lemma 14.71 we will rewrite the joint Laplace transform of ^ 1)1 Z{n 2 ,m 2 ) iii ^ form with 
structure closer to that of the Airy process (ITO . 

Proposition 4.8. 


= EE 


(- 1 ) 


m+n 


—' n!77i!(27ri)'"+"' 
n=0 m—0 ^ ^ 


'>0 “^>0 
n m 


A{x,t-,w,v) B{x,t-,w,w) 

C{x,t;v,v) D{x,t-,w,v) 


TT Trjvk - Wk) TT 7 r(v^ - Wi) -p-r -p-|- r(l + 7 - - Wfc)r(l + 7 - V£ - t’fc) 

sin(7r(ufe - Wk)) sin(7r(v^ - w^)) 1^1 r(l + 7 - - Ufe ) r(l + 7 - - Wk) 

=: V V /W 

/ V / ^ ^7n,n 


X 

k 

712 TUI 


n—0 m—0 


where in the above contour integrals the v and v variables run over Cs^ while the w and w variables 
run over £5 and the block matrices in the (n + m) x (n + m) determinant are given by 


A(a:,x;w,v)y = 


X w. 


/r( 7 -u;,)\™^ fT{vj) 


" Vr(7-7); Vrw 


1 < i, j < n. 


wi-^ wj- j r(7 - r(7 - Wj)”i 

X U 9 ^ Ml ^ ^^-^- — - 




1 < i < n, 1 < j < m. 


C(x,t;v,v),j = 


^ ^ r(7 — r(7 — ' 


L»(a;,T;w,v)„- = 


X u^ 


r( 7 -w,)V^ ^r(v,) 


r(7-v*)y V^K) 


1 < i < m, 1 < j < n. 


1 < i, j < m. 
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Proof. We begin with the formula obtained in Proposition l4.5l and use the identity zT^z) = r(l-|-z) 
to write the product 


nn 

fe=i i=i 


r(7 -wj- wfc)r(7 -vi- Vk) 
r(7 -Wi-Vk) r(7 - v£ - Wk) 


- TT TT (7 - - Tfc) (7 - - Wk) A A r(i -I- 7 - - wfc)r(i + 7 - 

(7 - - Wfe)(7 - - Vk) r(l -I- 7 - - tife ) r(l -I- 7 - - Wk)' 

Moreover, we write the term 

1 A TT{Vk - Wk) 


n 


sin(7r(pfc - Wk)) 


n 


n 


k=i ~ fe=i 


= (-1)” n / n 




fe=i 


7r(-Cfc - Wk) 
sin(7r(z;fe - Wk))'' 


and we also write the analogous formula for Ofci'^/“ '^l))^ by replacing in the above 
formula Wk by v^, and Tk by . Then we use Lemma 14.71 to obtain a block Cauchy determinant, 
to which we distribute into its first n rows the terms 
, r(7- 

" nwk)^- ’ 

into the first n columns the terms 


for /c = 1 ,..., n. 


{'Vk i)„Tk(vk-i) 


r(7 - Vk)""^ ’ 


into rows n + 1 to n + to the terms 


, (''< 2 ) „Te+n(ve-~f/2) b(vr) 


r (7 - Vr)"i ’ 

and into columns n + 1 to n + to the terms 

r(w^)"*i ’ 


-rt+n 


for fc = 1 ,..., n. 


for f = 1,..., TO, 


for ^ = 1,..., TO. 


□ 


We will further manipulate the formula obtained in Proposition 021 in order to bring it in a form 
suitable for asymptotics. Recall the definition of the functions G, F from (14.21) and denote by 

G(z) = G{z) - G( 7 / 2 ) and F{z) = F{z) - F{-f/2). 

In the block determinant that appears in Proposition 021 we multiply each column between 1 and n 
by the number exp (t 2 F('j/ 2 )) and each row between n +1 and to + r by the number exp (tiF{ j/2)). 
Moreover, we multiply each column between n -|- 1 and to + n by the number exp ( — tiF{j/2)) and 
each row between 1 and n by the number exp {^—t 2 F(^/ 2 )). This operation does not change the value 
of the determinant because if these number are pulled outside the determinant by mulitlinearity, 
they produce a pre factor exp {{nt 2 + mti)F{'^/2) — {nt 2 + mti)F{'^/2)') , which is obviously equal 
to one. On the other hand it provides a useful centering of the entries. Using the particular choice 
for values of rti, M 2 as in (02D and for the lattice points (toi, ni), (to 2 , 112 ) as in (14.11) and making a 
straightforward calculation employing Fubini, we have for fixed to, n that we can arrive at formula 

1^1 = (4.16) 


(- 1 ) 


m+n 


i\{2'!Tl) 

n 

xn 


m+n 


dr 


da; / dvdw / dvdw det 




A(a;,T;w,v) B{x,t;w,w) \ 
G(a;,x;v,v) D{x,t;w,\/) J 

r(i + 7 - - w/;)r(i + 7 - - Mfc) 

sin(7r(Mfe - Wk)) sin(7r(vfe - w^)) 4^4 1^1 ^ - Wk) 


7r(Mfe - Wk) 


/r"+'’ 


n 


77+ - Wi) 


nn 
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where 

A(a;,x; w, v)y = exp - G{wi)) +t 2 N'^^^{F{wi) - F{vj)) + r 2 N^^^{vj - Wi) + 

+ + f) (I - W*) - [Xj + Tj) Q - Vj'^ ^ , 

for 1 < *,J < ri, 

B{x,T;w,\N)tj = - exp ^ - N[G{wj) + G{wi)) + t 2 N'^^^F{wi) + tiN‘^^^F{v\ij) 

+ ( - Xj+n + T,) (^ - "'"i) + ( - ^j+n + Tj+„) Q - Wj) ^ , 

for i < i < n,l < j < m, 


G(a;,T;v,v)y = exp ^N{G{vj) + Giy^)) - t 2 N‘^^^F{vj) - F{vi) 

+ r2N^/^{v,-^)+nN^/^{s^,-l) 

— {xj + Tj) — Vj'j — (xj + Ti+n) ^ , 

for 1 < * < TO, 1 < j < n, 

D{x,T]w,v)^j = exp ^Af(G(vi) - G{wj)) + tiN^^^{F{wj) - F{vi)) + riiV^/^(vj - Wj) + 

+ {Xj+n + Xj+n) (^— — — (^Xj+n + Ti+n) (^— — Vi'j ^ , 

for 1 < i,j < TO. 


The use of Fubini is useful, in order to avoid crossing the poles at + Wfe =7 in the contour 
deformation that follows. Using the asymptotics (I2.15F we can justify the Fubini for fixed to, n (or 
for values that are a sufficiently small multiple of since the integrand in (14.1611 is bounded 

by 

iw^+^fei + ^ 

l<k<n l<^<m l<k<n 

l<£<7n 

_+n^ ^ Kl), 

l<i<m l<k<n 

which is integrable since m 2 — n 2 , ni — toi ~ 

The main asymptotic behaviour of each fixed summand in (14.1611 is captured by the determinant. 
The study of the asymptotic behaviour of the terms A, B, G, D in there is done via the steepest 
decent method. The basic idea of this method is to find the critical points for the integrand and then 
deform the contour of integration so that it passes through or close to the critical point and away 
from the critical point the integrands decay fast. The contours over which the v, w,v,w variables 
originally integrate are depicted in Figure [HI These contours need to be deformed to the steepest 
descent contours, which are depicted in Figure [TOl The contour G^, is symmetric over the real axis. 
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e 


7/2 


is 



Figure 10. Steepest decent contours 

Its upper part consists of a ray from 7/2 + , with fixed, small 17 > 0, leaving at angle tt/S 

until the point /2 + riN~^/^ and from this point it goes vertically to infinity. The Cy contour 

is also chosen to be symmetric over the real axis. It departs from 7/2 at an angle 27r/3, as a straight 
line until the point 76 ^’^'"/^, where it follows the circular arc ye''*, 27r/3 <t< 47 r /3 and, finally, from 
there it returns to 7/2 via a straight line. In order to motivate the choice of the contours Cy, Cy,, 
we need to identify the behaviour of the functions G,F around the (critical) point 7/2 and away 
from it. We do this by employing a Taylor expansion around 7 / 2 . To this end, we need to record 

G'{z) = ^'(z) + 4'(7 - z) + /.y , 

G"(z) = 'h'(z)-vh'( 7 - 2 ) , 

G"'(z) = ^"(z) + ^"(7 - z) . 


By the choice of = —24'(7/2) we have that G'( 7 / 2 ) = G”{'y/2) = 0, while G'"( 7 / 2 ) 7 ^ 0. 
Consequently, around the critical point 7/2 the function G(z) behaves as 

G(z) = G(z) - G( 7 / 2 ) = _ ^/ 2)3 + o{{z - ^/2f) (4.18) 

Moreover, 7/2 is also the critical point of the function F{z) because 

F'(7/2) = vI/( 7/2) - vIr(7/2) = 0, 

and around it we have 

F(z) = F{z) - F(7/2) = ^"y^^ (z - 7 / 2 )" + o((z - 7 / 2 )^), (4.19) 

The cubic behaviour of the function G as in (14.181) shows that along a straight line departing at an 
angle 7 r /3 from 7 / 2 , G is strictly negative with the value zero at 7 / 2 . Similar estimates hold along 
straight lines departing at an angle 27r/3. 

The behaviour away from the critical point 7/2 is described by the following lemmas. The first 
one is due to |BCR13j Lemmas 2.4, 2.5, which we recall : 

Lemma 4.9. There exists 7* > 0 such that for all 7 < 7* the following two facts hold: 
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(1) There exists a constant ci > 0 such that for all v along the straight line segments of Cy 

K[G(^;) - G(7/2)] < 5R[-ci73(z; - 7/2)3] ^ 

(2) There exists a constant C 2 > 0 such that for all w along the contour Cy, at the distance less 
than 7 from 7/2 


?ft[Giw) - G(7/2)] > 5 R[-c27^(w - • 

(3) There exists c > 0, such that for all v along the circular part of the contour Cy 


K[G(z;) - G(7/2)] < -c . 


The cubic behaviour in (I4.18|l suggests rescaling around 7/2 by change of variables 


Vi = iv3/3(u^ - 7 / 2 ) , 

Wi = - 7 / 2 ), 


V, =1V1/3(v,-7/2) , 

Wi = - 7 / 2 ), 

(4.20) 

Xi = N~'^/^Xi , 

n = 



Proposition 4.10. Recall the term Im,n in (I4.1bl) . Setting the parameters (toi, ni), (m 2 , n 2 ) as m 
su and the parameters U ‘2 as in ^depending all on IV/, then as fNl tends to inflnitg 

converges for any fixed m, n to 

w, w) \ 
;w,v) J ’ 
(4.21) 


. ri. — 


(- 1 ) 


n+m 


n!m! (27ri)'"+” 


dr / dvdw / dvdw 


da; det 


A*(a:,T;w,v) 

G*(x,t;v,v) 


D*{x, T 


where the variables v and v integrate along the contour U and w and w 

integrate along the contour y fQj- any horizontal shift 77 > 0. The 

block matrix in is given by 


A* {x, t; w, v)y = exp ^ 


+ r2{Vj - Wi) - {Xj + Ti)Wi + (Xj + Tj)Vj^ , 


- r2Wi - riWj - {-Xj+n + Ti)Wi - {-Xj+n + Tj+„)Wj^ , 

^ exp (^(SJ +7) - - *.3^7 

+ r2Vj + riVi + (xj + Tj)vj + {xj + ri+„)vi^ , 

ir{x,T-,^,'V» = exp - *3) + -7) 

+ ri(Vi - Wj) - {Xj+n + Tj+ri)^j + {Xj+n + T)+„)Vi^ . 


Proof. Let us start with the formula for 


(N) 


provided by (14.161) and denote 


CT 


m,n 


n m 


=nn 


r(i + 7 - - Wfc)r(i + 7 - - tifc) 
r(l + 7 - - rfe ) r(l + 7 - - Wk) 


(4.22) 


Step 1. (Contour deformation) We deform the contour Cs^ on which v and v are integrated, 
to contour Cy , as well as the contour on which w and w are integrated to contour Cw. Some care is 
required when deforming the w,w integrals. Let us look at it. For fixed values of r and x, the only 
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poles of the integrand (as a function of v, w, v, w) are those due to the cross term CTm,n coming 
from the poles of the gamma functions r(l -|- 7 — — Wfe) and r(l -|- 7 — — Wfe) when 

^ 'T 

j - We - Wk = {-- Wf) + {-- Wk) =-r - 1, for r = 0 ,- 1 ,- 2 ,... 

and 

1 Vk = (^-ve) + (^ - Vk) =-r - 1, for r = 0 ,- 1 ,- 2 ,... 

The deformation from to Cy crosses no poles, as on both contours Ji (7 — ve — Vk) > 0. 

Similarly, the deformation from is to Cw does not cross any poles, since on or on the left of the 
Cy, contour it holds that 5 R( 7 —w^ —uife) > —7 > —1, for 7 sufficiently small. Finally, the deformation 
is justified by suitable decay bounds near infinity between the initial and final contours. 

Step 2. (Localization around the critical point) Thanks to Lemmaand estimate (14.171) . 
the contribution to the contour integral away (e.g. outside a ball of radius 7 from the critical point 
7 / 2 ) is exponentially negligible as N tends to infinity. Then, again thanks to Lemma [4.91 we can 
further localize the contour integrals in a window of size around the critical point with a 

small cost e, uniformly in N. We now do a change of variables as in (|4.20|) to rescale this window 
by By using the Taylor expansion around critical point (|4.18|) and (|4.19ll . we can easily show 

that the contribution inside this ball is uniformly bounded by a constant, which will allow to pass 
in the limit in the quantities below. As TV —>■ 00 , we have (for fixed m, n) that 

” ™ r(l-7V-i/3(w, + ^i^))r(l-7V-i/3(v,+Ffc)) 

Mil r(l - N-y^{we + VU )) r(l - iV-i/3(v, + w^)) ^ ’ 


and moreover, 


ttN ^/^{vk - Wk) 
sin (7r7V-i/3(u^, - w^)) 


and 


■kN ^/^(vfc - Wfc) 
sin {jeN-^/^iyk - Wfe)) 


Using the Taylor expansions (14.181) and (I4.19|) for G and F, we have that in this limit and with the 
change of variables (14.201) . A,B,C,D converge to A*,B*,C*,D*, respectively, while the diagonal 
parts of Cv,Cy, are replaced by contours and {e“’^''/^IR>o+(5}U{e+’^''/^lR>o + 

(5}, respectively. □ 


By the multilinearity of the determinant, it follows that Im,n can be written as 

(^^ ('’’*+"> U )) m X n (B+n,Tj+n)) 


n!TO!(27ri)'"+"' i 


dx det 


+« 

>0 


(4.23) 


where 


II 

dvdic 

/ dx A*{x,T;w,v)ij, 

1 [R>o 


^ dfcdw 

/ dx B*{x,T;w,w)ij, 

J lR>o 

C (T^-j-yj , Tj ) — ^ 

^ dudv^ 

1 dx C*{x,T;v,v)ij, 
IR>o 

^ ("^z+n ;'^j+n ) — ^ 

^ dvdw 

/ dx £)*('?, f;w,v)y . 

' R>o 
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Let US, now, simplify the block determinant by doing the following change of variables 
2 


y/3 F"(7/2) 

G'"(7/2)V "G'"(7/2)’ 


V i-A - 


( 


( 2 n1/3„ F"h/2) 

W I—>■ (-; j w — t2- 

V G'"h 2 )J 


1/3 F"{^/2) 

V — 6l - 


W !->• — 


G'"(7/2)/ ^G"'(7/2)’ 

/ 2 n1/3_ F''h/2) 


G"'(t/2)^ ■ ''G"'(t/2)’ V G"'h/2)/ ‘Gh/2)’ 

Define = (—C 2 = — ^ ~ — c"'%^/ 2 ) ■ ^ ^“i® fairly straight¬ 

forward computation using the above change of variables and the integral representation of the Airy 
function (I4.4I1 . we obtain that 

pOC 

A'{Ti,rj) = / dxAi{clr 2 c^tl x n) Ai{clr 2 -\-c^tl x tj) , 

Jo 

poo 

B'{Ti,Tj^n) = - / - X + r*)Az(c7r2 + - X + Tj+n) , 

poo 

C'{Ti^n,rj) = / Ai(cJri-\-C2tl-\-x-\-ri^n)Ai(cJr2-\-C2tl-\-x-\-Tj) , 

^0 

poo 

D'in^n^rj^n) = / dxAi{clri-\-c^tl-\-x-\-ri^n)Ai{clri-\-c2tl-\-x-\-Tj^n) • 

The only tricky part of the computation that leads to the above expression is to take care of some 
multiplicative factors, which are independent of the variables w, v, w, v, S and appear in the entries 
of the blocks A'{n, tj), B'{n, Tj+n), C'(Ti+n, tj), D'{n+n, Tj+n)- Denoting these terms by exp(Gj), 
it turns out that Cij are equal to 

, . r ,2 (J^"(7/2))^ ^ f"(7/2) , _ f"(7/2) , 

^{i<2<"} 1 (1^2 2 I ^ *G"'(7/2) ^ 


3G'"(7/2) J G"'(7/2) 


, ^2X"(7/2))^^ ^^"(7/2), _F"(7/2)^ 

^{n+l< 2 <n + m} ] I lc\\ I 


3G"'(7/2) ; G"'(7/2) 


'G"'(7/2) 


. , ^2X"(7/2))7 ^^"(7/2), F"{^/2)^ 

l{i<i<n} S — I f’2 — *2 ) '7v7i7Z7TA\^^ ~ D'^^7777 — ~'t 2 


3G"'(7/2) j G'"(7/2) 


G"'(7/2) 


“1“ ^{n+l<j<n+m} 


/ ^, (F-(7/2))^ \ ^^^(7/2) . _ i^"(7/2) . \ 

V ^ ^ 3G'"(7/2) y G"'(7/2) ^ ^G'"(7/2) 7' 


These terms can be pulled out of the block determinant, by multi linearity, and they then cancel 
out. 

An immediate change of variables brings us to expression (14.91) . thus bringing us to the end of 
the road map. 
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